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Abstract. This article concerns arbitrary finite heteroclinic networks in any phase space dimen- 
sion whose vertices can be a random mixture of equilibria and periodic orbits. In addition, tangencies 
in the intersection of un/stable manifolds arc allowed. The main result is a reduction to algebraic 
equations of the problem to find all solutions that are close to the heteroclinic network for all time, 
and their parameter values. A leading order expansion is given in terms of the time spent near 
vertices and, if applicable, the location on the non-trivial tangent directions. The only difference 
between a periodic orbit and an equilibrium is that the time parameter is discrete for a periodic 
orbit. The essential assumptions are hyperbolicity of the vertices and transversality of parameters. 
Using the result, conjugacy to shift dynamics for a generic homoclinic orbit to a periodic orbit is 
proven. Finally, equilibrium-to-periodic orbit heteroclinic cycles of various types arc considered. 

1. Introduction. Heteroclinic networks in ordinary differential equations or- 
ganise the nearby dynamics in phase space for closeby parameters. They thus act as 
organising centres and explain qualitative properties of solutions, and predict varia- 
tions upon parameter changes. This makes heteroclinic networks a valuable object 
in studies of models for applications. When all vertices in the network are equilib- 
ria much about such bifurcations is known. Recently, heteroclinic networks whose 
vertices can also be periodic orbits have found increasing attention. 

This article concerns the unfolding of hnite heteroclinic networks consisting of 
hyperbolic equilibria and periodic orbits in an ordinary differential equation 

-^u(x) = f(u(x);fx), (1.1) 

with x <G R, u(x) € M." and parameter /x £ M. d for arbitrary n and sufficiently large d. 

Any solution that remains close to the heteroclinic network of (jl.lj) assumed at 
/i = for all time can be cast in terms of its itinerary in the heteroclinic network. 
See Figure ITTTI for a simple example. For any itinerary, bifurcation equations will be 
derived for the locus of parameters of all corresponding solutions to The idea to 

formulate the unfolding in this way is borrowed from previous studies of heteroclinic 
chains of equilibria [551 EH] • 

Bifurcation studies from heteroclinic chains with equilibria mainly concerned ho- 
moclinic orbits and generated a huge amount of literature, see, e.g., [51 ITU1 UM [TBI I2o1 
[3H [371 EI] j us t to name a few, and heteroclinic loops between two equilibria, see, e.g., 
[30I3[12[T1[I3[371[4T1[42]. Heteroclinic cycles with periodic orbits have found 
increasing attention recently, see, e.g., j^lH[r^|^|^|g^|5^l55]. 

The main new contributions of the present work are rigorous results allowing for 
periodic orbits in general heteroclinic networks, for tangent heteroclinic connections, 
and to formulate the bifurcation equations in a general form that can be used as the 
basic building block for a specific study. It is hoped that this makes the results useful 
for readers with applications to specific cases in mind. 

A large field of applications are travelling waves in evolutionary partial differential 
equations in one space dimension whose profiles solve an equation of the form (|l.ll) , 
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Fig. 1.1. Sketch of an itinerary for a 2-homoclinic orbit, (a) The homoclinic orbit q-[ (solid) 
and its asymptotic state p* at jX = 0, and a 2-homoclinic orbit (dashed). The Poincare- section 
£ is just for orientation here, (b) Schematic plot of the itinerary (solid) for a 2-homoclinic orbit 
(dashed) where 92 = 93 = <?i> Pi = P2 = P*- 



but also for instance Laser models arc often reduced to this form. There is a very 
large amount of such analytic and numerical studies in the literature, e.g., [5I lll?Il2"2l 
US [3U [351 [361 [381 [39] to hint at some. 

In the applied literature such bifurcation equations arc frequently derived formally 
by a geometric decomposition in terms of local and global maps, e.g., [U [TU 151] . 
The justification in particular of the local map is an issue and, if linear, requires 
non-resonance conditions on eigenvalues or else dimension dependent normal form 
computations. Other issues are the form of parameter dependence and the persistence 
of solutions upon inclusion of the higher order terms of the original vector field. These 
problems do not arise in the approach taken here, and the results can provide a 
rigourous foundation of formal reductions. 

The reduction to bifurcation equations in this paper is motivated by the so-called 
'Lin- method' described in [55], which is a Lyapunov-Schmidt reduction for boundary 
value problems of the itinerary. This method has been used and modified in a number 
of ways and contexts for equilibria, e.g., [T71 [T|)J 2U] . Tangent intersection of stable and 
unstable manifolds have been considered mostly for homoclinic orbits, i.e., homoclinic 
tangencies, a paradigm of chaotic dynamics; Lin's method in this context has been 
used in |16j . Periodic orbits introduce technical complications and for Lin's method 
these have been overcome in |32] and, using Poincare- maps, in |33| . Transversality 
studies with respect to parameters in related cases were done in [13]. An ergodic 
theory point of view is taken in [I] [TT1 [531 (Ml HZ1 HH] > and further papers by these 
authors, looking for instance at properties of non-wandering sets. More recently [3J 
treated periodic orbits in a very promising way using Fcnichel-coordinates. 

Here 32 is used as a starting point, and equilibria or periodic orbits as vertices 
are treated in an essentially unified manner. Symmetries or conserved quantities are 
not used, but a generic setting is assumed. In contrast to [32], winding numbers of 
heteroclinic sets are not considered, and the underlying heteroclinic network is held 
fixed. Together with [32] this exposition is self-contained, but somewhat technical, 
and parts of [32) have to be repeated and improved in order to track higher order terms 
in this extended setting. The precise statement of the main result Theorem 14.31 can 
only be given rather late after a number of preparatory steps, notation and definitions. 
In particular, this includes £|3J where we obtain suitable coordinates near the vertices. 
We next describe the main result, and refer to $5] for sample applications. 

1.1. Description of the main result. For a chosen itinerary the method is 
a Lyapunov-Schmidt reduction which yields algebraic equations that relate system 
parameters fj, to certain geometric characteristics Lj , Vj at each heteroclinic connection 
qj that the solution follows, perhaps repeatedly, and which connects vertex Pj~\ to 



Pj. The time spent between Poincare sections at qj_i(0), and Ej at <7j(0) is 2Lj 
for Lj 6 [L*, oo ) if is an equilibrium. If pj is a periodic orbit, this time is in general 
only approximately 2Lj since we normalise Lj <G {lTj/2 :/6N, £Tj/2 > L*} for the 
minimal period Tj of vertex j. For tangent hetcroclinic connections or more than one- 
dimensional hcteroclinic sets, un/stable manifolds of Pj-i and pj have a more than 
one-dimensional common tangent space at (fj(0), and Vj are the coordinates on that 
space, except the flow direction. The location of Lj , Vj in the itinerary is illustrated 
in Figure 11.21 

The system parameters fi g R d can generically be assumed to unfold each hete- 
roclinic connection by a separate set of parameters /x* £ K. dj which, however, must 
coincide at repeated connections. Here dj is the codimension of the jth connection 
and d = dj without repeating the same connection in the sum. The geometric 
characteristics couple these parameter sets, but to leading order only to the nearest 
neighbours j ± 1. If dj = for all j, then all hetcroclinic connections are transverse, 
and the result proves the existence of solutions for any itinerary, and an expansion 
for the coordinates in the Poincare-sections. Otherwise, an expansion of fi* in terms 
of v r , L r is provided as described next. 

Let Kj S and aj be the real and imaginary parts of the leading un/stable eigen- 
value or Floquet exponent at vertex j. For 7 e R r let Cos(7) = (cos(7i), . . . , cos(7 r )). 
For all j with dj > 1 and for sufficiently large min r { L r } and small sup r { | v r \ } , there ex- 
ist Pj , 7j € K dj ' , and linear maps (3j , 7^ , /?" , 7" , Q , , (j , £J- , (" , , as well as quadratic 
maps Tj (zero if tj =0), so that a solution follows the chosen itinerary, if and only if 



$ = TjM + e-^^aJLj + ^ivj^Lj+^qivj^Lj^) 
+e -^U L i^Cos(2a S j_ 1 Lj_ 1 +j;(vj- 1 ,Lj_ 2 )q(vj- 1 ,Lj_ 2 )+n j 



Here /i* = whenever the j and j' in the itinerary correspond to the same actual 
heteroclinic connection. The coupling to the nearest neighbours is given by 

%(v,L) =/3j+/3'jV + p f jBf +1 (L), 

7 *(v,L) =i j+7 'jV+^B*j_ 2 (L), 

Q(v,L) =Cj+CjV + q'Bf +1 (L), (1.3) 

Q(v,L) = &+£» + #BS_ a (L), 

B S r /U (L) = C- 2 < /UL Cos(2^ /U i + fi /U )Cr /U , 

where (3 s r /u G R dr and Cr /u arc linear maps. 

Negative Floquet multipliers, i.e., negative eigenvalues of the period map, have 
Floquet exponent with imaginary part ft/Tj. Since L = lTj/2, £ £ N, the argument 
in the cosine terms is tt£ which generates an oscillating sign as £ is incremented. 

Note that if the itinerary has repetitions, then Vj,Lj have to satisfy solvability 
conditions from repeating the corresponding equations. Each repetition yields new 
parameters Vj,Lj, but all other quantities in (( 1 . 2|) are the same if the underlying 
heteroclinic connections is the same. 

The significance of (|1.2[) lies in the order of the remainder term IZj, which, for 
certain <5"^ s > 0, and arbitrary S > 0, is given by 

\vj\ 3 + e' 2 ^ L J (e- 5 i L * + \v j+1 \(e^- K ^ L i+ i + \v j+1 \) + e (^-^ +1 )L J+1 ^ 
+c -2k=_ 1 l j _ 1 ^g-a^z^ + | Uj ._ 1 |( e (<5-«5-2)ii-2 + \vj_ 1 \) + e {S - 3KB i-2'> L i- 2 ^ . 
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Fig. 1.2. Schematic illustration of an itinerary near index j, and the location of [i*,Vj,Lj. 
The arrows indicate the flow direction, connections can be copies of the same actual heteroclinic 
connection. 



In particular, IZj is higher order with respect to at least one of the cosine terms if 
Lj ~ L r> r = j - 1, j - 2,j + 1. 

The application of the above result to a Shil'nikov-homoclinic orbit to an equi- 
librium yields the same bifurcation equations as in |26j . and as in that paper, many 
of the seminal results by Shil'nikov [37] follow from leading order analyses. Note that 
the resonant be treated by the above result as well. See [6] for 

resonance at homoclinic bifurcations. 

Remark 1.1. Under the 'flip' condition Rank(C 3 ) = or Rank(£j)j = for 
the leading eigenvalue, the next leading terms need to be taken into account for an 
unfolding. Viewing Q , £j as parameters, the bifurcation of solutions can be understood 
from TheoremU^if 2k] < 2««._ 1 + for Q = or 2«S_ 1 < 2k] + Sj for ^ = 0. 
To overcome the barriers involving Pj_i and p" to the next order eigenvalues requires 
a more refined setup beyond the scope of this article. See \3J$ and also \2b\ H2J/ 
for such considerations in case of equilibria applied to homoclinic bifurcations. 

The main result unifies the treatment of equilibria and periodic orbits as ver- 
tices of a network: for the reduced equations the only difference between equilibria 
and periodic orbits is that Lj is a semi-infinite interval for an equilibrium, but the 
above defined discrete infinite sequence for a periodic orbit. The discrete sequence 
essentially counts the number of rotations that the solution makes about the periodic 
orbit. Note that replacing a periodic orbit by an equilibrium has consequences for the 
codimensions of heteroclinic connections. 

The reduced equations for a specific case can be determined in three steps. First, 
choose the itinerary of the solution type of interest. Second, determine the codimen- 
sions, including tangencies, of all visited heteroclinic connections. Third, copy the 
equations from Theorem 14.31 for each element in the itinerary with positive codimen- 
sion, and remove geometric characteristics that do not occur according to the type of 
itinerary and tangencies. In case of repetition in the itinerary, the locus of parameter 
values for the solutions should be found by analysing the arising algebraic solvability 
conditions (which can be highly non-trivial). Similarly, in case of tangencies, the locus 
of turning points or folds can be determined. 

To illustrate this and the applicability of the abstract results, some sample ap- 
plications for specific heteroclinic networks are presented in $5] In particular, for a 
generic homoclinic orbit to a periodic orbit conjugacy to (suspended) shift dynamics is 
proven. In addition, equilibrium-to-pcriodic orbit heteroclinic cycles of various types 
are considered, and 2-homoclinic orbits are studied for the first time in this context. 

Note that the main result separately concerns each solution type as encoded in 
each itinerary. This is suitable, for instance, when looking for the aforementioned 
travelling waves. In some cases a whole class of solutions or even the entire invari- 
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ant set can be characterised directly. However, our results do not provide stability 
information of the bifurcating solutions or hypcrbolicity of invariant sets, or other 
ergodic properties. Sec [35] for stability results in homoclinic bifurcations using Lin's 
method (where additionally PDE spectra are considered). As mentioned, the above 
result does not provide an expansion in general in the case of vanishing leading order 
terms ('flip conditions'). 

This paper is organised as follows. Section [2] contains details about the setting 
and some preparatory results. In [J3] a suitable coordinate system is established for 
trajectories that pass near an equilibrium or periodic orbit or lie in un/stable man- 
ifolds. The main result is formulated and proven in $4] Finally, £[5j contains sample 
bifurcation analyses and illustrate how to use the main result. 

Acknowledgement. This research has been supported in part by NWO cluster 
NDNS+. The author thanks Bjorn Sandstede, Ale Jan Homburg and Alan Champ- 
neys for helpful discussions. 

2. Setting and Preparation. The basic assumption is that at fi = (jl.ljl 

possesses a finite heteroclinic network C* = (C^C^) with vertices p* £ C* being 
equilibria or periodic orbits, i £ I\, and edges q* £ C%, i £ 1%, being heteroclinic 
connections. Rather than unfolding C as a whole, we consider the following paths 
within C separately. Here we set J° := J \ min J (where min J = if there is no 
minimum). 

Definition 2.1. A possibly infinite set C = ((Pj)j£ji(qj)jeJ°) with pj £ C*, 
j £ J and qj £ C\ , j £ J° , is called itinerary if for all j £ J° the edge qj is 
a heteroclinic connection from Pj-i to pj (or homoclinic if Pj-i = Pj) and either 
J = { — oo, ... 0, 1}, J = {1, 2, . . . oo}, J = Z or J = {1, . . . , m} for some m > 1. 

For ease of notation we say that a sequence yj of 'objects ' ( numbers, vectors or 
maps given in the context) with j £ J° has reducible indexing (with respect to C) if 
Vj = Vj' whenever qj = qy for j,f £ J°. 

Note that an itinerary can cycle arbitrarily and perhaps infinitely long within the 
heteroclinic network viewed as a directed graph, and the labelling can differ from that 
in C* . Any itinerary has a (possibly non-unique) reduced index set J re( j C J so that 
qj / qj< as well as pj ^ p r for j ± j', with j,f £ J r ° cd := J rod n J°, and j,j' £ J rod , 
respectively. Associated to this is C rcc i = ({Pj)j£j Tt , d , d ) C C, which may not 

itself be an itinerary (though it contains one) . 

Let Jg C J be the index set of all equilibria pj in C and Jp = J \ Je that of all 
periodic orbits. We set J£ := J E H J° and Jp := J P n J°. Finally, Tj > denotes the 
minimal period of pj for j £ Jp, and we set Tj = for j £ Je- 

In the following, unless noted otherwise, we consider an arbitrary fixed itinerary 
C. However, until S]only neighbours qj and qj-i are relevant. 

For j £ J let &j(x, 0) = Aj(x)c FjX be the Floquet representation of the evolution 
of the linearization v — d u f(pj(x);0)v of (|l.ip in pj. Here v = dv/dx, and the 
matrices Aj(x) satisfy Aj(0) = Id, Aj(x + Tj) = Aj(x) for j £ Jp and x £ R, and 
A(x) = Id for j £ Je (in which case Fj = d u f(pj\ 0)). 

The basic assumption about (|1.1[) and the heteroclinic network is 
Hypothesis 1. The vector field f in (|1.1[) is of class C k+2 for k > 1 in u and fx. 
The equilibria or periodic orbits p* , i £ I\, are hyperbolic at fi = 0, i.e., for any C the 
matrices Fj have no eigenvalues on the imaginary axis, except for a simple eigenvalue 
at the origin ( modulo 2iri ) if j £ Jp . 
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Here a simple eigenvalue has algebraic and geometric multiplicity one. Hypothe- 
sis [T] implies that the spectrum spec(Fj) of Fj has the un/stable gaps 

ft" := min{Re(i/) : v G spec(Fj), Re(u) > 0} > 
k S j := -max{Re(^) : v G spec(Fj), Re(i>) < 0} > 0. 

Since C\ is finite, the gaps Kj are uniformly bounded from below in j G J. For 
convenience we choose arbitrary Kj > 0, j € J with reducible indexing, such that 
Kj < min{«;^,ft"}. We also need the gap to the next leading eigenvalues/Floquet 
exponents. Let is r , r = 1, . . . , n be the eigenvalues of Fj and define 

Pj := min{|Re(iv)| — ftj : Re(v r ) < 0, Re(i/ r ) ^ Kj, r = 1 . . . , n}, 
Pj := min{Re(i/ r ) — k" : Re(^v) > 0, Re(V r ) ^ ft", r = 1 . . . , n}. 

Leading stable eigenvalues of a matrix are those with the largest strictly negative 
real part, and leading unstable those with the smallest strictly positive real part. For 
the main result, we will assume that these leading eigenvalues are simple as expressed 
in the following hypotheses. 

Hypothesis 2. Consider the leading stable eigenvalues or Floquet exponents 
at pj. Assume that this is either a simple real eigenvalue Vj or a simple complex 
conjugate pair Uj, Vj with Im(i^) ^ 0. 

Hypothesis 3. Consider the leading unstable eigenvalues or Floquet exponents 
at pj. Assume that this is either a simple real eigenvalue Vj or a simple complex 
conjugate pair Vj, Vj with Im(^j) ^ 0. 

To emphasise where these hypotheses enter we will not assume them globally, 
which has the effect that a priori exponential rates for estimates are not k^ u , but 

ft^/ u — Sj for an arbitrary Sj > due to possible secular growth. In the following 
Sj denotes a priori an arbitrarily small positive number, which may vanish under 
Hypothesis 03 El 

Hence, for suitable Xj G [0,7}) as well as asymptotic phases ctj G [0,Tj) of qj 
with respect to Pj we obtain the estimates (see, e.g., [7]) 

\qj(x)-pj(x + aj)\ < Ce&- K V x , x > (2.1) 
\q j+1 (x + Xj) -pj(x + ctj)\ < Ce (5 ^^ )|x| , x < 0, 

where C > depends only on qj,qj+\ and Sj. For j G Jp the requirement (|2.1[) of 
equal asymptotic phase for qj(0) and qj+i(xj) determines Xj up to multiples of Tj 
and uniquely in [0, Tj). For j G Je we have Pj{x) = pj and set ctj — Xj =0. 

To distinguish in- and outflow at pj we denote q~j(x) := qj + i(x + Xj) for j-1 S J°. 

Let &j(x, y) denote the evolution of v = d u f(qj{x);0)v and &j(x,y) that of 
v = d u f(qj(x);0)v. Hypcrbolicity of pj gives the following exponential dichotomies 
for j G and trichotomies for j G Jp (see, e.g., [32]). 

Notation. Indices separated by one or more 'slashes' as in indicate al- 
ternative choices for the statement with all these indices chosen equal at a time. 

There exist projections ^ a J c ^(x), continuous in x > 0, and ^^ c ^ u (.t), continuous 
in x < 0, such that the following holds. Set ^°^ u (x,y) := $j ;(x, y)^ c ^ u (y) and 
<& S j (x, y) := &j(x, y)^ S j^ C ^ U (y), respectively. 
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. For J G Jp : Rg(^(x)) = span 

• For j G Jg : = yj/c = 0, 

• The projections arc complementary: + \pu + = id, vfi s (vl/ u + = 0, 
* u (* s + = 0, * c (* s + = 0; analogous for #® /c/u , 

• The spaces Rg( x I ,s (x)) and Rg(^ u (a;)) are unique, the spaces Rg( v I /U (a:)) and 
Rg(4' s (x)) arbitrary complements such that the previous holds, 

• The projections commute with the linear evolution: 

$f /u (x,y) = ^f /u (x)^(x,y) and l»f /u (x,y) = ¥/ c/u (x)^(x,y), 

• They distinguish un/stable and center direction: there is C > depending 
on 5j > and qj , <jj such that for all u G M. n 



(2.2) 



® s j(x,y)u 


< Ce {s i 


-<®\*-v\\u\ 


a; 


>y>0, 


$j(x,y)u 


< Ce {5 i 


-*?)l*-wl| u | 


• y 


> x > 0, 


® c j(x,y)u 


< C\u\ , 


x,y > 0, 






$j(x,y)u 


< Ce (5j 




, X 


<y<o, 


&j(x,y)u 


< Ce ( ^ 


-«5)k-y|| u | 


u 


< x < 0, 


®j(x,y)u 


< C\u\ , 


< 0. 







We denote the un/stable and center spaces, respectively, by 

E«' s ' c (x) :=Rg(*f ^'/ s / c (x) :=Rg(*f / c (z)). 

Definition 2.2. _Fbr a decomposition E® F = M" we denote by Proj(E,F) 
M. n — > WL n the unique projection with kernel E and image F. 

In order to link the trichotomies of the in- and outflow near pj , we define 

P;(L) := Proj(£»(Z) + Ej(-L), E*(-L)) 

Pf{L) :=Proj(E» c (-L),EJ(L)) 

P?(L) := Proj(^(i) + Ej(-L), E](-L)). 

We also define the aforementioned sets of travel time parameters 
Kj(L) : 



[L, oo ) for j G Je 

{£Tj/2 : iTj/2 >L,^GN} for j G J P . 



It is shown in [32], Lemma 2, that there is L° > such that for L G Kj{L°) 
the P^ c ^ n (L) arc complementary projections, = for j G J£, and the norms 
are uniform in L. 

In order to control the leading order terms in the bifurcations equations we make 
the following change of coordinates locally near all pj, j G J ro d- In the new 'straight' 
coordinates the un/stable manifolds locally coincide with the un/stable eigenspaces 
of the linearization in pj, respectively. For periodic orbits the strong un/stable fibers 
locally coincide with the un/stable eigenspaces. Since these are graphs over the 
eigenspaces and tangent at the equilibrium or periodic orbit this is straightforward. 
See e.g., (3.27) in [33]. However, as in [33] this change of coordinates is an obstacle 
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to apply the method within the class of semilinear parabolic partial differential equa- 
tions. However, in [5] this problem has been circumvented in a way that should apply 
here as well. 

To emphasise the effect of this coordinate change and to make the notation of 
estimates throughout the text more readable, we define for j <G J and any Sj > 0, and 
5j = if explicitly mentioned, the terms 



e ( s j- K j) L 3 m straight coordinates, 

&. s i~ K i) L i a priori, 
c ( s j- K j) L i m straight coordinates, 

and set Rj = R 3 ■ = for j $ J. 

Notation. In the following we use the usual order notation a = 0(b) if there 
is a constant C > such that \a\ < C\b\ for all large or small enough b and norms 
as given in the context. In terms of Lj this is always as Lj — > oo. In a chain 
of inequalities for such order computations we allow the constant C to absorb other 
constant factors and take maximum values of several constants without giving explicit 
notice. 




The next lemma is the basis for estimating error terms in the following sections. 
Lemma 2.3. There exist C > and L 1 > L° depending only on qj,qj and Sj 
such that for LL 3 G K{L r ) the following holds for all j G J° . 
1. \P™{L 3 )^{L 3 )\ < CR 3 , \Pj C (Lj)^J(—Lj)\ < CR 3 . 
f \P™(L 3 )(p 3 ia 3 + L 3 ) - q 3 (L 3 ))\ < CR 3 R 3 , 
■ \\Pf(L 3 ){p 3 {a 3 - Lj) - q 3 (-L 3 ))\ < CR 3 R 3 . 

3. The above holds under Hypothesis with 5 3 — in R 3 and under Hypothesis 
[5| with Sj = in Rj . 

4- Under Hypothesis [H or [J respectively, there are vectors v^ s in the leading 
un/stable eigenspaces of Fj such that 

v S j ^ <^> lhnswpliL(qj(x))/x — —tfj, 
vj ^ <^> limsupln(gj(a;))/x = Kj, 

IE— > — oo 

and such that under Hypothesis® and for any 6j < mhi{K B j,Pj} we have 

$5(0, -L)P?(L)( P j(aj +L)- q 3 (L)) = c 2F > L v s 3 + 0( e - 2 ("S+*5> £ ), 

and under Hypothesis^ and for any SJ < min{K", p"} we have 

^(0,L)P?(L)(pj{ aj — L) — q 3 {-L)) = e~ 2F ^v] + 0(e~ 2 ^ +s ^ L ). 

Proof. For readability, we set ctj = 0, see (|2.ip . Let ^j' u (x) be the stable/center 
and unstable eigen- or trichotomy projections on the whole real line, x £ K, of 



v = d u f(pj-0)v, 
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which trivially exist by the Floquct form. Note that these di/trichotomics differ from 
those of the linearization in q 3 . 

1. First note that 

*§(£) =Proj(£;| u (L),^(L)) 
P™(L) = Proj(%(-L), ££(L) + EJ(—L)), 

so that PJ u (I/)^-(L) is determined by Ej(L) — Ej(—L). (An appropriate norm 
for estimating this difference goes via suitable bases of these linear spaces.) Since 
for L G Kj{L°) it holds that **•(—£) = for all j G J° (the projections arc 

constant for j G Jg) we have 

E*(L) E*(-L) = E*{L) Rg(^(i)) + Rg(tf°-(-Z,)) - E*(-L), 

and we will estimate the two differences on the right hand side separately. 
General perturbation estimates of dichotomies (e.g. Lemma 1.2(i) in [34]) imply 

- &j(L) = 0(e^- K ^ L ) and - = 0(e^~ K ^ L ). 

Hence, on the one hand E^L) - Rg(^.(L)) = e>(e (5 ^ )L ). 
On the other hand we can write 

*5(-£) - &j(-L) = Id - *™(-L) - (Id - *f (-£)) = *™(-i) - * 3 cu (-i), 

and, due to asymptotic phase we have ^(?j(— L) — -§^Pj{— L) = 0(e^ Si ~ K i^ L ) in 
the center direction. Therefore, ^(L) - ^(-L) = 0{c {S] -^ ]L ) and Rg(^(L)) - 
E)(-L) = 0(e^- K ^ L ). (And analogously - ^(L) = 0(e^~ K ^ L )). 

In combination, since Ej(—L) C kerP"(L) the weak version of the first estimate 
follows. The strong version of this estimate in straight coordinates is a consequence 
of the fact that then E*(L) = Rg(**-(L)) for all L > L 1 for sufficiently large 
L 1 > L . Hence, for L > L\ we have 

E*(L) - E](-L) = Rg(^(-i)) - E](-L), 

which implies the stronger estimate. 

The proof of the second estimate is completely analogous. 

2. Since the stable manifold is a (at least) quadratic graph over the stable eigenspace 
for j G Je and center-stable trichotomy space for j G Jp at p 3 we have that 

Pj (L) - q 3 (L) = i>f(L)(p 3 (L) - q 3 (L)) + 0{{p 3 {L) - q 3 {L)f). 

On the other hand, as in the proof of the previous item, we can replace ^ S 3 (L) by 

&j{L) with error of order e( 5j ~ K ^ L . Since E?(L) lies in the kernel of Pj n (L) and 

p 3 (L) — q 3 {L) = 0(e( Si ~ n i' L ) there are non-negative constants C* and C* such 
that 

P™(L)( Pj (L) q 3 {L)) = P™(L)^(L)(p 3 (L) - q 3 {L)) (2.3) 

+0(C,e^-"S-' , ?) £ + C*c 2 ^- K ^ L ). 

For j G Jp recall that the asymptotic phase as L — > oo of q 3 (L) is that of Pj{L), 
hence q 3 {L) lies in the strong stable fiber with phase L. Strong stable fibers are 
(at least) quadratic graphs over the (strong) stable trichotomy spaces so that 

& 3 (L)(p 3 (L) q 3 {L)) = 0((p 3 (L) q 3 {L)f) = 0(c^~^ L ). 
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The claimed weak estimate follows from combining this estimate with (|2.3|) . 
The stronger estimate for straight coordinates means C* — in (|2.3|) . Indeed, 
since the strong un/stablc fibers and (center) un/stable manifold coincide with 
their tangent spaces at pj the higher order corrections disappear and (|2.3I) holds 
with C* = 0. 

The estimate for Pj(oij — L) — qj{—L) is completely analogous. 

3. For simple eigenvalues or Floquet exponents there is no secular growth so that the 
rates of convergence are in fact the leading rates. 

4. This is a reformulation of results in [33] as follows. 

Concerning the right hand sides without <&^(0, — L) and $"(0,L), respectively, 
Lemma 10 in [32] yields the expansions 

P°(L)( Pj ( aj + L) - qj (L)) = ML)^^ + 0(e- 2 K+^), 
PfiLypjiaj — L) — fc(-L)) = Aji-L^-^v] + 0{e- 2 W+W L ), 

for certain u"^ s in the leading un/stable eigenspace of Fj. Note that in the present 
case a and v from that Lemma are constant, and L G Kj(L ) so that Pj{oij +L) = 
Pj(aj — L) and Aj(L) = Aj(—L). 

Concerning %(0, — L) and $"(0,L), Equation (5.7) in [33] shows that for any v 
there are w u / s in the un/stable eigenspace of Fj such that 

*5(0, -L)v = c F ' L A J {-L)- 1 v s + O U-W +8 P L } , 
*j(0,L)v = e-^AjiL)- 1 ^ + O (e-^H"^) . 

Since Aj(L) — Aj(—L) for L 6 KjlL 1 ), and these expansions only depend on the 
error to the asymptotic vector fields, combination of this with the previous step 
proves the claim. 

□ 

3. Coordinates of trajectories. Following and improving [32], in this section 
we establish a suitable coordinate system for the (n— l)-dimcnsional set of trajectories 
that pass nearby pj. We consider the difference V = (w,w) between solutions u to 
and <7j, qj, where u(0) ~ 9j(0) and u(2L) ~ 9j(0). (In this section 9j can 
be any orbits that lie in the stable and unstable manifolds of pj, respectively.) We 
determine all such V by an implicit function theorem, where the time shift along a 
trajectory is removed to make V unique. For equilibria this is done in the usual way 
of [26] by imposing certain boundary value data of u in Poincare sections attached to 
the in- and outflowing solutions qj and <jj, and adding a continuous parameter L so 
that 2L is the time spent between the section. 

For periodic orbits it would be natural to do the same, only L would not come 
from a connected unbounded interval in order for V to be small. However, due to our 
approach via a certain variation-of-constants solution operator, we slightly deviate 
from this, see also [35]. Briefly, in order to control the integrals over the centre part 
of the trichotomy projections, we use exponentially weighted spaces. This causes 
difficulty to control the centre projection of another integral term from coupling in- 
and outflow, which stems from the deviation of the phase with respect to pj at x = L. 
Due to asymptotic phase this term vanishes in the limit L — > oo, but integrability 
requires a good estimate. We avoid this and at the same time remove the phase shift by 
simply requiring that the centre parts of w{x) and w(x) at time x — L, x = — L vanish. 
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This precisely disallows time shifts and the result is equivalent to the approach by 
Poincarc sections. In particular, the resulting V give a parameterisation of all orbits 
near pj with normalised travel time sequence. A posteriori, the reconstructed solutions 
approximately satisfy boundary conditions in suitable linear Poincarc sections and 2L 
is the approximate travel time between these. 

Notably, this parametrises trajectories in a neighborhood of {qj(x) : x > 0} U 
{(jj(x) : x < 0}. An alternative to this approach is to parametrise trajectories in a 
small neighborhood of pj and then insert a 'global' trajectory piece between the inflow 
and outflow boundaries of these neighborhoods, see [21 HI 157]. 

3.1. Passage coordinates. We choose the boundary value data in subsets of 
Poincare sections defined via the trichotomies. For j G Jp the space (<7j(0), qj(0)) + 
Ej(0) x -^"(O) is (n — l)-dimensional since the flow direction counts towards stable 
and unstable directions for periodic orbits. On the other hand, for j G Jg we have 
qj(x) G Ej(x) and q(x) G Ej(x) so that after removing the flow directions for in- and 
outflow n — 2 dimensions are left. As mentioned, this is compensated by a continuous 
'travel time' parameter L G K^L 1 ). To eliminate the flow direction for j G Jg we 
define 

Q) ■= Proj (Ej(0) + span(q,(0)),/^(0) D spanfe (0))^) , 
Q] ■■= Proj (E}(0) + span(4(0)), EJ(0) D span(^ (0)) x ) . 

Here E 1 - is the orthogonal complement of a linear space E. For the outflow we define 
Qj accordingly, and set Qj ■= Q] ■= Id for j G Jp. 

Now boundary data in Q^E^O) x Q^E^(0) for all j G J° excludes precisely 

the flow directions at qj(0) and qj(0). Recall that Ej(0) and £^(0) are determined 
uniquely by the di-/trichotomy. 

Finally, for all j G J we define the Poincare sections 

Sj :=^(0) + Q5£5(0)+^(0), 
tj :=^(0)+£5(0) + Q^(0). 

Hence, for j G we solve the boundary value problem (jl.lj) subject to u(0) G 
Sj and u(2L) G Sj. As mentioned, for j G Jp there are technical reasons not to 
consider this boundary value problem. In addition, the set of travel times has to be 
'phase coherent' in order for the variations V to be small near the periodic orbit. 
For convenience we a priori restrict L to the set Kj(L ) which already appeared in 
Lemma l2~3l For general Poincare sections, the set Kj would need to be adjusted and 
therefore, in general, differs for each n, which is inconvenient for the leading order 
expansion of parameters. The tradeoff is that orbits starting in Sj do not necessarily 
lie exactly in Sj for L G Kj^L 1 ): instead of having zero center part at time 2L we will 
require this at time L. Since near qj(0) the flow acts as a diffcomorphism between 
any choice of hypcrplanes transverse to the flow this is no restriction, but rather a 
normalisation of the discrete travel times. 

Due to the geometric interpretation, we refer to the boundary data as coordinate 
parameters and denote 

:= Q^ s (o) x </;!■:;;■ ui. 
ii 



Fig. 3.1. Illustration of the notation for the flow near a periodic orbit pj with minimal period 
Tj , j G Jp. Theorem \3.2\ see also Corollarv \3.3\ show that trajectories near {qj{x) : x > 0}U {q«(x) : 
x < 0} are parametrised by G E|(0), G E^(0) and L G {iTj/2 : I G > L }. 

with elements uj = for G Q^(0), u)j G Q]EJ(0), see FigurcO 

Since we also want to parametrise un/stablc manifolds including the flow direction 
for equilibria, we define in addition 

fij :=E*(0) x £^(0). 

We now turn to the aforementioned difference V between solutions and hetero- 
clinic orbits. Given a solution u for parameters \x we call any 

V(x, x; fi, u, a, L) := (u(a + x) — qj (x),u(a + x + 2L) — qj(x)) 

a j-variation of u and denote the components as V = (w,w). A solution u of () 1 . 1 [) 
for a; € [0, 2L] can be reconstructed from a given j-variation if on the one hand 

^w(x) = d u f(qj{x);0)w(x) + gj(w(x),x; fi), (3.1) 
-^w(i) = d u f(qj(x)\0)w(x) + g 3 (w{x),x; fi), 

for a; G [0, L] and x G [— Z, 0], where 

gj(w{x),x;n) ■= f(Qj(x) + w(x);fi) - f(qj(x);0) - d u f{q j {x)- 1 0)w(x), 
gj(w(x),x;fj,) := f \qj{x) + w(x); fi) - 0) - d u f(qj(x);0)w{x). 

On the other hand, the reconstructed orbit is given by 

q 3 {x) + Wj(x) , x G [0, L] 

& (x - 2L) + ibj (x-2L) , x G [£, 2i] 

and is continuous, if w 3 (L) — Wj(—L) = —qj(L) + q~j(~L). Therefore, we define 

:= qj(—L) — qj(L), 

which will be the main contribution to the expansion given in (|1.2[) . 

We look for solutions that are simultaneously close to qj and qj, so that the 
j-variation V is small and given by an implicit function theorem. Since bj{L) is 
asymptotically periodic for j G Jp as L — > oo the aforementioned 'phase coherence' 
condition appears, and for L G Kj{L) we indeed have smallness: 

6j(L) = Pi ( aj + L) - qj (L) + qj{-L) - p 3 -(ay + L) = 0(c^ L ). 
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Moreover, Lemma T2.3f 4) implies for Lj G Kj^L 1 ) that 

Pf(L j )b j {L j ) = 0{R j ) (3.2) 
Pf{L i )b j {L j ) = 0{R j ). (3.3) 

The trichotomies imply for x, — x G [0, L] that 

R" - E s {x) x £ c (a:) x £ u (a:) - E s (x) x i? c (z) x E u (x) 

and provide a decomposition w = w *■ + ufi + w" into Wj (x) := ^ (x)w(x) and 

u) = w S j +w^ + w^ into w a J c ^{x) := $^ c ^ u (x)t&(:c). We use the analogous superscripts 
for the decomposition of gj and gj . 

Similar to [34] , Lemma 3.4, the following estimates hold. We first change coordi- 
nates and rescale time so that (see [32]) there is £o > such that for \fi\ < e we have 
Pj(x; 0) = Pj{x; fx), j G J re d, i.e., pj is independent of \i for \p\ < eq. 

Lemma 3.1. There is C > depending only on Sj and qj,qj such that for 
x, — x > 0, 

\ gj (w,xw)\<c(\w\ 2 + \M\w\+eto-W*)) (3-4) 
\$(w,x;im)\ < C ((1^1 +e^-"^ I )H a + W\w)\ + (3.5) 
\gj(w,x;n)\ < C (\w\ 2 + \fi\(\w\ + e^~ K ^)) (3.6) 
\gj(w,x;v)\ <c((\w»\+ c^^\)\w\ 2 + | M |(|^| + e^"^ 1 * 1 ) (3.7) 

Proof. By definition of straight coordinates, for w G R", we have ^>j(x)f(pj(x) + 
^(x)w\0) = 0, where ^ (x) arc the projections of the linear evolution at pj as 
in the proof of Lemma [231 From that proof we know ^-(x) — &j(x) = 0(e^ 5i ~ K ^ x ); 
recall \Dy (x) are the projection of the trichotomies $(x, y), see (|2.2[) . Therefore, 

^(x)f(p 3 (x) + ^(x)w,0) = 0(e&-*}>). 
In the following we omit the argument x for readability and set / s := tyj(x)f. Hence, 

d U uf s ( Pl + w;0) = 0(|u$| + e^-^ x ). 

We have 



9j{w, x] n) = f{q 3 +w;n)- f(qf,0) - d u f(q 3 ]0) 



w 



f(qj + w;ii) - f{q 3 + w; 0) + /(<& + W, 0) - /(<&; 0) - d u f{qf,Q)w 

9 M /(q + w; r/i)/idr + / / d uu f(q + tsw; 0)rw 2 drds, 
o Jo Jo 



and d[_if(pj +w\fi) = C(|w|) since f(pj]fi) = f(pj;0). It follows that d^f = 0(\w\ + 
e ( s j- K j) x ^ which proves the first claimed estimate. We also infer d^f s = 0(\w s \ + 
C ( s j- K j) x ^ which proves the second claimed estimate. The proof of the remaining 
estimates is completely analogous. □ 
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Based on this, V can be found with uniform estimates in L for j £ Jp in the 
weighted space 

X rhL = (C°([0,L],R") x C°([-L, 0],M"), || • \\ n , L ) , 



{W,W)\\ ri . L 



II' 



\\w\\+ L = S up{\e*i*w(x)\:xe[0,L]}, 
\\w\\- L =su V {\e^w(x)\:xe[-L,0]}. 

By Lemma \3. 11 for any < r)j < there is a constant C > independent of L 
such that 



<C((\\w\\+ >L ) 



ImI), 



||ff i (ti>,-;M.i)ll^ l£ <c((IHI~,,l) 2 + Im 



(3.8) 



The coordinates of trajectories will be those (uij,(j,,L) with L £ Kj(Li), ujj = 
(uipUJ^) £ fij that generate a fixed point of the map 

g j (w,w;uJ :j ,n,L) : X VtL -> X^z,, j £ J°, 

defined next. The maps Qj can be derived from a variation-of-constants solution of 
(|3.1[) decomposed suitably by the trichotomies and are given by 



Gj(w,w\Uj,iJ,,L) 



( *5MK \ 

+ Jl ®T( x ' y)9j(w(y), y; n)dy 

+^(x, L)Pf{L){c,{L)uj + w; ^, L) + 

$y(sc,0)w/ 

+ Jo ®](x,y)9j(w(y),y, v)dy 

\ -$-(i,-i)/- c (L)(c J (i) W + c),Ku;; Ai , J L) + fo J (L)) / 

where the horizontal line separates first and second component Qj = (Gj,i, Qj,2)- The 
terms coupling these components are 

tj{L) Uj = l>, u (-J,o)^ - -i-^ /..();.,;. 

p — L pL 

dj(w,w;^,L) = / $](-L,y)gj(w(y),y;(j,)dy- y)gj(w(y), y; n)dy. 

Jo Jo 

By Lemma 12.31 there is C > depending only on 5j and qj , qj such that 

\P?(L)cj(L)uj\ < C(Rje^- K ° )L \^\+e^- K * )L \u>]\) (3.9) 

(3.10) 



P S j{L)tj{L)uj\ < C(e^- K ^ L \L0 S j\ + R jC ^- K ^ L \uJ\) 



Note that ^(L)G jt i(x = L) = and ^(-L)£j, 2 (i = -£) = 0. As mentioned 
above, for j £ Jp, this is equivalent to fixing the phase on a reconstructed orbit from 
a fixed point of Qj. It is shown in |32j . Lemma 4, that fixed points of Qj indeed 
generate the aforementioned reconstructed orbits u(x) of for a; £ [0,2L]. The 

following theorem proves that all orbits can be obtained in this way. Recall that 
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ojj G flj contains the flow direction for j G Je so that the map from fixed points to 
trajectories at x = x = is not injective, while for ujj G fij it is. 

Notation. B(X,p) := {x G X : \x\ x < p} where ■ \ x is the norm of X. 

The following theorem provides the coordinates of trajectories near qj, cjj (in fact 
near any in- outflow pair at pj); this is formulated more explicitly in Corollary 13.31 
below. 

Theorem 3.2. Assume Hypothesis^ and take j G J, as well as any rjj G (0, Kj) 
if j G Jp, and rjj G [0, Kj) if j G Je- There exist e > 0, L* > L 1 depending only on 
qj,qj such that the following hold for L G Kj(L*), p G B(M. d ,e), LOj G B(Clj,e). 

1. There exists a unique Vj — Vj(u)j, p, L) G X^.l such that Vj = Qj(Vj\Wj, p, L). In 
addition, Vj is C k smooth in (uij,p) and for j G J£ also C k smooth in (ujj,p,L). 

2. Let V(x, x; p, u, 0, L) = (w(x), w(x)) be the j -variation of a solution u of (|f .![) with 
u(0) G Tjj and, for j G J^, u(2L) G Sj. If \V(0; p, u, a, L)\ < e then there is unique 
a* = 0(\Pj(L)w(L; p,u, a)\) such that 

V(-;p,u,a*,L) = Vj{w h p,L) 

foruj = (Q5*5(0) W (0;^u,a*,i),g;'^(0)^(0;^ U ,a*,L)). 

3. For any 8j G (i]j,Kj), there is C > such that a fixed point (W, W){uij, p, L) of 
Qj(-,ujj, p, L) satisfies, for x,—x G [0,i], 

||W(^,/i, J L)||+ >i <C(| M | + |u; J s |+^), 

||W(^-,/i,i)||- >i <C(| M | + |^|+4.) ) 
\^ s J {x)W J {uj J ,p,L){x)\ KCeV'-W*, (3.11) 
\^(x)Wj{cjj,p,L){x)\ <Ceto- K W. (3.12) 

4- Under Hypothesis^ or\3[ the estimates (|3 . 1 1 1) and (|3.12p hold with 5j = rjj in Rj 
or in Rj, respectively. 
Concerning the required smoothness of /, the loss of two degrees of differentia- 
bility is due to the coordinate changes (which can be performed simultaneously) that 
involve /' and that g contains /'. 

Before proving this theorem, to emphasise the coordinate system character we 
reformulate part of Theorem 13.21 in the following corollary taking ojj G flj . Recall 
from the beginning of this section that then the set of parameters (u>j,L) is (n — 1)- 
dimensional for all j G J°. The theorem in particular shows that this is also true for 
the set of fixed points. 

Corollary 3.3. For any j G J° there is e > and L* > L 1 as well as a 
neighbourhood U of \qj(x) : x G [0, L]} U {(jj{x) : x G [-L, 0]} such that the following 
holds. The set of solutions u of (jl.ip that lie in IA and satisfy u(0) G Ej, and, 
if j G Je, u(2L) G Sj for L G Kj{L*) is in one-to-one correspondence with the 
parameters {(oJj, p, L) : u>j G B(Qj,e), L G Kj(L*), \p\ < e} of fixed points of Qj. 

Proof of Theorem \3.2\ Items 1 and 2 are a consequence of Theorem 1 in [32] as 
follows. There it was assumed that \p\ < ee _,,ji , but in the present case we can take 
\p\ < e due to the following estimate. For any < rjj < Kj we have 

\\d w9j (w,-,p)\\+ L < tfidH+i + H), 
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and the same with hats in the || • ||~ L -norm (see [33] Lemma 3.1). This estimate 
improves the corresponding estimate at the end of the proof of Lemma 5 in [32] as 
needed. All constants are uniform in L £ Kj{L*). 



Items 3 and 4 improve the estimate of Theorem 1 in [32] using Lemma 12.31 as 
follows. Theorem 1 in [32] states 

\\(w 3 ,w 3 )\\, hL < c(H + H + c^ L ), 

where for j G Je we can set 77 = since there is no centre direction. 

To improve this we consider Wj and Wj separately and note that the coupling 
of these only enters through Dj, while the dependence of Wj on and Wj on uj* 
is only by Cj. On account of (|3.9p , (|3.10p . the decomposition of ui in the claimed 
separate estimates of Wj and Wj follows. The estimates (|3.2[) . (|3.3[) prove the claimed 
separation for the remainder term bj(L). Again note that 77 = is allowed for j 6 Je- 

It remains to estimate Oj and to prove the pointwise estimates for W| = QjW 
and Wj 1 = <& S jW. We first consider the pointwise estimates. By definition of Gj,i 



W s 3 {x) = *5(s,0)«5 + / ^(x, 2/ ) 53 (I^(y), 2 /; M )dy, 
Jo 

and so using (|2.2[) and (|3.5p . for a weight 77 to be determined there is C > such that, 

e vx \W s 3 (x)\ < C L^- K >\^\ + J* e^- K > +K jV ((\W?(y)\ + e^- K V y )\Wj(y)\ 2 

+ \v\(\W 3 (y)\+e^-^y))dy) 

< C fe (lJ -"S> !B |a>$| + jf e ^- K i )x+K i* ((e~ w || WJ|| + L + e ( ^"^ )!/ ) 

x e- 2 ™(||^||,t L ) 2 + |Ml(e- w ||W7l|+ £ + e^-"5>»)) Ay) . 

Taking maxima over x this implies for any < 77 < K S j (and < 77 < Kj for j £ Je) 
that 

||^||+ L < c (k s | + (1 + llWJII+JdlWfll+J 2 + , 

and in particular for any 5 > there exists C > such that 

II^IIJ-^l < C & \W](x)\ < Ce^*, x e [0,L]. 

The estimate for Wj(x) is completely analogous, only now Kj bounds 77. 

We now turn to the required estimate involving j. Substituting p. lip . (|3.12p 
into p.5p and (|3.7p gives C > such that 

\g S j(w,x W )\ <c(e^-^ x (\w\ 2 + M)) (3.13) 

|#(<MiAOI <c(e^-«?)l*l|(|i2>| a + |A i |)) . (3.14) 
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Concerning these estimates yield (suppressing y,fi in gj) 



-L 



| / &(-L,y)gj(Wj)dy\ < C / e-W^e^-^dW^y)] 2 + \p\)dy 



-L 



^ L e- 5 ^\\\W j \\l L f6.y + |m|) 



< C 



(e-^dl^ll-^ + H 



Analogously, 



L (ll^lltL) 2 + H 



Hence, using Lemma |2"T 31 there is C > depending on 8j > such that 



Im|), (3-15) 
Im!)- (3-16) 
□ 



IP/Ci)^,^;^)!^ ^(11^11-^ + ^(11^11+^ 

|Pf(L)^-(^,^;L,/i)| <C(RjRj(\\Wj\\-. tL ) 2 +Rj(\\Wj\\+ tL ) 2 

This completes the proof of the claimed estimates. 

Note that the differences in this result between an equilibrium (j £ Je) and a 
periodic orbit (j £ Jp) are twofold: 

f . The ranges of L values measuring the time spent near pj are a semi-infinite interval 

for j £ Jg and a discrete infinite (phase coherent) sequence for j £ Jp, 
2. In the periodic case the exponential weight rjj must be strictly positive. 

3.2. Stable and unstable manifolds. For homoclinic or hctcroclinic connec- 
tions it is helpful to parametrise the stable and unstable manifolds of pj in the same 
way as above. This is in fact simpler than the general case and we can set L in the 
definition of Qj to infinity, which gives, for j £ J°, 



*|(«,0)w5 

+ SL ®T( x > y)9j ( w (y)' v\ ^ d y 



( 4>](x,0)cjJ \ 

+ J^oo y)g~j(w(y), v, v)dy 
\ + So ^']( S: ^y)9j(w(y),y, fi)dy J 



The same change in the proof of Theorem 13.21 gives the following corollary concerning 
the parametrization of un/stable manifolds W s / U (pj). 

Corollary 3.4. There exists e > such that for all ({j,, Uj) £ M d x fi*. with + 

< e the operators Q°° and Gj° have unique C k smooth fixed points W°° (x; tu^ , fi) 
and Wj°(x;u)j, fi) satisfying the following. 

1. Forj 6 Je the maps lo* Wj°(0; Wj, /i)(0) andu:^ i— > Wj°(0; U)J, fi)(Q) parametrise 
W s (pj) andW n (pj) near qj (0) and qj(0) over Ej(0) andEJ(0), respectively. 

2. For j £ J P , a £ [0,Tj) the maps uj] i-> W°°{a; lu s v fi) and luJ m- W?°(a;uif, fi) 
parametrise the strong stable and unstable fibers of pj with phase ctj +a over Ej(0) 

and EJ(0), respectively. 

3. Theorem\?[]^3) holds for W°° with Rj = and W°° with R 3 = 0. 
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Fig. 4.1. Schematic illustration of adjacent j -variations. 



4. Bifurcation equations. Based on the results of the previous section we 
derive reduced equations whose solutions parametrise all solutions of (jl.ip that are 
near the chosen itinerary C. Throughout this section we take uij = (wj, u/J) G fij. 

In order to reconstruct solutions of (JTTTJ from the variations (Wj,Wj) about 
adjacent qj these need to fit together continuously. By definition of the variations 
this means solving (up to the flow direction as shown below) the system of equations 

Wj(x = 0;u$,&?,n,Lj) = Wj-i(x = 0;w5_ ls ^_ ls (4.1) 

where j G J°. System (|4.1[) is closed if J = Z, and closing conditions are required for 
J with upper or lower bound. 

In case of finite J reconstructed solutions are either heteroclinic from p\ to p m 
('het' in short) and, for p% = p m , homoclinic to p\ ('horn') or periodic orbits ('per'). 
Note that the same periodic orbit is a solution for any periodically prolonged itinerary, 
and in this case we implicitly assume C is a heteroclinic cycle. The remaining cases 
are semi-unbounded J for which we require the corresponding solution to lie in the 
un/stable manifold of pj with the largest or smallest index, respectively. 

More formally, this means 

'het': ©(O) + W 2 (0) G W u ( Pl ) and <? m _i(0) + W m -!(0) G W s ( Pm ), i.e., 

Wa(0;wa,^L a ) = T^r(ai;w?,/«), 
W m _i(0;w m _i,/i,i m _i) = W^(a m ;uf m ,ij), 

'horn': (really the same as 'het', here extra only for clarity) 

<Z 2 (0) + W 3 (0) G W u ( Pl ) and g m _i(0) + W m ^(0) G W s ( Pl ), i.e., 

^ m _i(0;o; m _i,/i,L m _i) = Wi°(ai;u>l,[i), 

'per': Wi(0; ux } M> L i) = W m (0; u m , A 4 ; L m ), 
'semi-': q (0) + W o (0) G W s ( Pl ), i.e., 

^o(0;w ,^Lo) = Wr(ai;ul,ti), 

'semi+': ^(0) + W 2 (0) G W u ( Pl ), i.e., 

^ 2 (0; W2 ,M,i 2 ) = VK 1 oo (a 1 ;^ , , M ). 

In order to unify notation for these cases, we set L\ = oo for 'semi±' and L\ = 
L m = oo for 'het' and 'horn' so that all equations are of the form (|4.1j) . We thus omit 
the superscript 'oo' and take indices modulo m + 1 for 'per'. System (|4.1[) then needs 
to be solved for j G J° with the modified definition 

!J mod m + 1 for 'per', 
J\ {1} for 'semi-', 'horn' and 'het', 
J for 'semi+' and when J = Z, 
18 











U i+ i 


— > — 








— ^ 













9i(0) 



*+i(0) 



Fig. 4.2. Schematic plot of the flow box near the trajectory between <7j(0) and <jj_|_i(0). 
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Free travel time parameters are then Lj with j G J L where 

L J J° for 'semi±', 'per' and if J° = Z, 



{2, . . . , m - 1} for 'horn' and 'het'. 



The parameters of fixed points of (Qj)j^jo are thus LJj, j E J°, and Lj, j e J L , and 
the actual system parameters fi S In the Lyapunov-Schmidt reduction we first use 
the coordinate parameters ujj, and then, if needed, express the system parameters n 
through the time parameters Lj and possibly remaining coordinate parameters. This 
also determines the generic minimum number of parameters needed for the unfolding 
of the part of the network that is visited by the selected itinerary. 

If C contains a sequence of adjacent periodic orbits, the requirement of equal 
asymptotic phase in Theorem 13.21 for in- and outflow at each of these may require 
different qj + i(0) and <jj(0), i.e. Xj ^ in the definition of qj. In that case solving 
(|4.1[) requires a nontrivial shift in the flow direction. However, this direction is not 
directly available since we removed the flow direction from the coordinate parameters 

To trivialise the matching in this direction we change coordinates locally in a 
neighbourhood of the trajectory segments Yj := {qj(x) : < x < maxjTj, Tj_i}} 
for all j G Jp, to obtain 'flow box' coordinates so that the flow is parallel to Yj in a 
tube about it, see Figure Since C* is finite we can choose a uniform tube radius. 
Note that this change of coordinates is independent of the changes near pj performed 
above. 

Remark 4.1. 

1. Since fixed points of Qj are coordinates of trajectories ( Corollary \3.3\) there is a 
bijection between solutions (up to time shifts) of system (|4.1|) with all closing con- 
ditions and solutions of (jl . lj) that stay in a certain neighbourhood of the itinerary, 
if we require minimal period for periodic solutions. 

2. Due to the flow box coordinates near problematic qj(0), it is in fact not necessary 
to solve (|4.1[) in the flow direction qj(0): Even if the orbits reconstructed from fixed 
point components Wj of Qj and Wj-\ of Qj-i do not fit together in the flow direc- 
tion, a unique trajectory is selected, see also Figure \4-2\ Proof: By construction, 
all trajectory segments 



Uj(x) 



<&(aO+Wj(a:), xe[0,L,} 

qj {x- 2Lj ) + Wj(x-2Lj), x<E [L 3 ,2Lj], 



j G J L , are continuous at Lj. A priori we would require the jumps Uj(2L) — Uj+i(0) 
to vanish. However, since the vector field in the flow box is parallel to qj the 

19 



coordinates in the un/ stable trichotomy spaces Ej^ 1 1 (x) and Ej(x) do not change 
within the flow box. Therefore, the segments Uj anduj + \ lie on the same trajectory, 

if and only if their j- and (j + 1) -variations have same coordinates in (0) and 

E s / U (0). 

Recall the spaces of coordinate parameters if? := QjEj(0) and HJ := QjEJ(0), 
and that these do not contain the flow direction. 

To motivate the following definitions, notice that and explicitly appear 
in (|4.1[) when substituting the definition of Qj at x = x = 0. In fact, (|4.ip projected 
onto Ej := Hj + HJ_ l can be solved using uij, and we therefore define 

Pj : I'nOE: .E,). 

Lyapunov-Schmidt reduction now consists of solving the system (|4.1|) projected first 
by Pj and then substituting the result into the projection by Id — Pj . In this process 
the flow direction need not be considered as shown in Remark 14.11 Therefore, the 
directions that are unreachable by coordinate parameters are 

E):= (^(0)+i^(0))\ 

Hence, dj := dim(E^) is the number of reduced equations at ^(0) that need to be 
solved by system parameters. 

Definition 4.1. Let J b c J° be the set of indices for which dj > 1. 

We set d := X^eJt>nJ d ^ anc ^ snow that this is the number of parameters 
needed to unfold C rc d and thus to locate the solutions selected by the choice of C. We 
call additional parameters auxiliary. 

If Hj fl is non-trivial, then the representation of Rg-P,- by ufj + is 

not unique. To make it unique, we remove the intersection from Hj and denote the 

remaining coordinate parameters by Vj € Hj n H^_ 1 . More precisely, we define 

Pj := Proj([ff? n H^,H] n H^) (4.2) 
H] := ker(P,) n H), 

so that for ujj £ Hj there exist unique Vj £ Rg-Pj and Qj £ Hj with ujj = Vj + uij. 
Definition 4.2. Let J 4 C J° be the set of indices for which dim(Rg(Pj)) > 1. 

We denote the collection of all these coordinate parameters by 

and endow V with the sup-norm. Parameters Vj occur if the tangent spaces of stable 
and unstable manifolds coincide in more than just the flow direction. A transverse 
heteroclinic set of two or more dimensions occurs for j £ </' \ J b , which means that 
the 'linear' codimension 

n+l- dimW s (pj) - dimW u ( Pj --i) 

is negative and gives the generic dimension of tangency minus the flow direction. 
This only uses information from the un/stablc dimensions at the asymptotic states, 
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Fig. 4.3. Illustration of the notation for a tangent heteroclinic connection. 



and tangency of the manifolds may be higher dimensional, and can also occur for 
positive linear codimension. The above defined dj includes this by accounting for the 
intersection of Hj and HJ, and is always larger than or equal to the linear codimension. 
We therefore refer to dj as the codimension of qj. Note that transverse heteroclinic 
connections occur for j £ J° \ J b and tangent directions transverse to the flow (a 
tangent heteroclinic connection) occurs for j £ J' fl J b , see Figure l4~3l 

To capture the leading order effect of parameter variations on the j-variation we 
define the Melnikov-type integral maps for j 6 J b 

Mj : R d -> Ej , ft ^ / 0)a* , a jtr (y)) dy a% r , 

r=l ^ R 

where a° r £ It", r = 1, . . . , dj is a basis of E}j with reducible indexing, and a,j :r (y) is 
the solution to the adjoint linear equation 

a = -d u (f(q j (y);0)) t a, 

with a(0) = a" r . On account of (|2.2I) _M_, is well-defined. 
Note that auxiliary parameters ft lead to a modified map 

(ft, ft) i ^ Mf/x + Mj/i G -E b . 

The complete Melnikov-map for c? parameters is then 

M : R d -> £ b := ] [ £ b , // h> (M^) jeJb . (4.3) 

Hypothesis 4. ker(.M) = {0}. 

Repeated elements in C mean repeated rows in M.. For an equation M.ft = X this 
means that the solvability conditions on the coordinates Xj of X in i? b are Xj = Xj> 
whenever qj = gy, j,j' £ J b . To solve the remaining part of A4/i = X separately 
in each £" b as far as possible we change parameters as follows. Under Hypothesis 2] 
the Melnikov map M. Te d ■ ^ d —> I\je.J b n.J rad s invertible. Set ft = M lc dft and 
f(u;fi) := f{u;M.~ eA ft) and omit the 'check' in the following. In the new parameters 
A4 re d is the identity on E^° in the sense that 

dj 

r=l 

for a unique subcollection fij = {fij r ) r =i....,d j of parameters, and fi = {jlj)j£j Ted nJ b - 
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To unify notation of bifurcation equations for parameters and solvability condi- 
tions we define itinerary parameters fi* for all j G J h as follows. Set p* = fij for 
j G Jrcd H J b and. for j G J h \ J re di = Aj' whenever j' G J b is such that qj = qy . 
Due to the above change of parameters, solutions to Mjfi — X can be cast simply as 

/'.; -v,. ./ . •/''. 



In preparation of the main theorem statement we define for j £ J PI J* the map 

which measures the quadratic separation of the tangent manifolds by Vj G Rg(Pj) 
and 

On account of (|2.2[) Tj(vj) is well-defined and Tj(vj) = 0(\vj\ 2 ). The following terms 

will capture the leading order effect of the neighbouring itinerary elements and give 
rise to the expansion of the bifurcation equations. 

ir.L s ;: -I-;;!). /. ; i/>; ; ;/ v ^;/ v i. 
in/.,:: /,,:/>;; 

Prom (1221), (O and flH3]) we infer 

1^(^)1 - 0(eV- H V L * Rj) = 0(R 2 3 ), (4.4) 
|^(Lj)| = ©(e^-"^^^) = O(flj). (4.5) 

The following hypothesis concerns intersections of the heteroclinic orbits and the 
spaces Ej as well as Ej with leading un/stable fibers and trichotomy spaces, respec- 
tively, and excludes flip bifurcations. 

Hypothesis 5. Let lA be the leading stable /unstable eigenvalues or Floquet 



exponents at pj with Im(z/^ u ) G [0, 2ir), and k^ u := Rc(z/^ u ). 
ln(gj(x)) _ u 



ln(aj jr .u (x)) 



' 3 

1 . lim sup — V- * J v ;y = i% 

X^ — OG X 

2 . Umsup ln fa-iW) = ^_ x) 
5. 3r u G {1, . . . , dj} swc/i lim sup ■ 

X" — »• OG — X' 

^. 3r s G {1, . . . , dj} such that limsup ^ n ^ a ^ r ■. 

x—t — oo X 

ln($"(x,0)?;) 

5. 3v G Ej such that limsup = uj, 

2:— >oo X 

ln(^.(x,0)u) 

6. 3v G Ej such that limsup = 

x— > — oo X 
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To emphasise the local coupling in the itinerary and to show conjugacy to symbolic 
dynamical systems, sec H5.2.21 for an arbitrary < A < 1, we use the weighted norm 

||L|| C := £ A^l (4.6) 

J'GJ L 

on the space C(L*) := IlrGJ L Kr(L*) of the sequence of travel time parameters Lj. 

In the precise formulation of the main result given next, we identify with the 
isomorphic R d J and for j G J b the following denote linear maps, where M C R d J xd J 
are diagonal matrices. 

$ : Rg(P,'+i) -> R d ', 7- = Rg(^-i) -> R d ', 
: £^ +1 (0) -> R*i, i! : £f_ 2 (0) -> R dj ', 

0,fj 6 M, Cj : Rfi(P,-+i) M, : Rg(Pj-i) -> A/, 

g:^ +1 (0)^M,^:^_ 2 (0)^M. 
Here ker(/3/), ker(Cj') 3 kerP j+1 , ker( 7 j'), ker(£") D kerPjLi so p'jP J+1 = P'j etc. 
Let il = (1lj) je jb, and Cos(/3) = (cos(ft), . . . ,cos(/3 r )) for /3 G R r . 
Finally, u = (wjjjejo, Co = (Q s ,u") € := U, • //; x P^), and £ : n jeJ o £y, 
where (l and .E are endowed with the sup-norm. 

Recall that Lj = oo, if j g" J L and = if j £ J'. Similarly, we make the 
convention that a quantity vanishes if its label is outside its range. 

THEOREM 4.3. Under Hypotheses [7] and[7J for a given itinerary C with closing 
conditions, if required, there exist L* , £*,£* > depending only on C re d so that for all 
8 > the following hold. 

1. For all j G J h assume Hypothesis [H and [3| at p r for r = j, j ± 1, j — 2 with 
leading un/stable eigenvalues iv = —Kr + io> , respectively. Take any <5^_j < 
minf/c^, /o^jj and (5" < min-fre^, p"}. For r = j,j ±l,j - 2 set S r = 
if r E J h fl Je and otherwise 5 r = S. 

There exists unique Pj,^j G R dj and /3j, /?", 7^, 7", £j, ^, Cj> £j> Cj'; as above, 
as well as unique C k smooth w s , w u , ft) : B(V,£*)x£(L*) -> B(R rf x x £ b , e*) 
swc/i i/iai (/U, w s + v, Cj u ) solves (|4.ip /or j G J° if and only if (j,j(v,L) satisfy 
with p.3[) /or j G J b . quantities except (v,L) have reducible indexing and 



Hj = 0( \ Vj \ 3 + e- 2K "^ e-"^ + |%'+i|(P i+ i + |^+i|) + R 



-6"L 



_|_g — 2«j_ 1 ij_l 



— ft, = 0(e-( a *?+^)^ + e -( 2K 5-i+ 5 J-i)^-i + I jay 
duj 

Finally, Rank(Cj) > 1 under Hypotheses \9(1) and (3), Rank(£j) > 1 under Hy- 
potheses\B(2) and (4). The analogous statement holds for Cf+i> Cj—2- 
2. For j G J°, Sj = S and rjj as in Theorem \3.2l and R = J2ieJ b nJ d + ^» 

solutions to (|4.ip with \ (p,,uj)\ < e* and \L\ > L* satisfy 

m(v,L)\\l Lj - OQvjl + Rj + \v\ 2 + R 2 ), 
\\Wj(v,L)\\- jLj = OQojl + Rj + \v\ 2 + R 2 ). 

If j G Je, then under Hypotheses^ and [3] at pj we can take Sj — in Rj and Rj. 
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3. There exists a neighbourhood IA ofUj e j° d {qj(x) : x G R} such that the set of(L,v), 
\v\ < £», L G C(L*) for which there is a solution to (|4.1[) with \fi\ < £* is bijective 
to the following set of((J,,u) G B(R d , e*) x C°(R,R n ): u solves ([II]) withu(0) G £1 
and u{x) G IA for all i£K, and there exists (xj)j£j° C [0,T U ) with Xj+i — Xj > 
minimal such that u(xj) G Ej, where T u G R U {00} is the minimal period of u. 
In [J5] the use of this somewhat abstract result for concrete cases is illustrated by 

a number of examples. See §1.11 for a discussion of Theorem 14.31 



The remainder of this section is devoted to the proof of Theorem 14. 3[ which 
proceeds in the two Lyapunov- Schmidt reduction steps 1. solve (|4.1j) by the coordinate 
parameters uij, 2. solve the remaining equations except the flow direction by system 
parameters fx. 

4.1. Solvability by coordinate parameters. 

In Ej, the leading order dependence on ji will stem from 



M^-= / ^ 3 (y)d f ,f(q J {y);0)ndy, 
Jr 



which is well-defined due to (|2.2j) for $j 



Pj^(0,y) ,y>0, 
Pjty-^y) ,y<0. 



Lemma 4.4. There are C, L*, si, e%, £3 > depending only on C re d and Sj such 
that there exist unique C k smooth (uj s ,ui n ,1Z)(v,L, fx), (uj b ,w u ,1t) : B(V, £1) x C(L*) x 
B(WL d ,e 2 ) -> x E, that satisfy, for j G J°, 

1. P 3 (Wj(0; Q} + Vj,&j, n, Lj) - Wj-iiO; + Vj-^Qj^n, = 0, 

2. |c^ s | + |^_! - < C + + + -Rf-i), 

5". wj - = -Vj + M j fi- Bj(Lj) - B^^Lj-i) + Kj,!, where 

K jtl = O + i?; 3 _! + {Rj + Rj-! + \v 3 \ + \ft\)\ft\ + \vj\iRj-! + R, + \ Vj \)) 

4- For any solution of (|4.1[) multiplied by Pj that has < £2 and \ui\ < £3 there 
exists unique (v, L) such that u> — (uj b (v, L, fx) + v, uj n (v, L, jj,)). 
Note that Mj, -B", Bj_ x have reducible indexing. 

Proof. For j G J° we need to solve 

Pj (pjAW^Wj^^L^O) -6? i _ lj2 (W 7 -_ 1 ,^_ 1 ;^_ 1 ,M,^-- 1 )(0)) = 0. 

Let vj G Rg(Pj) and G for the decomposition of = Qj + Vj. Reordering 
terms, using Rg(Pj) = Ej, and that the center direction lies in the kernel of Pj this 
equation becomes 



u>* - = -vj + U 3 . 2 - Bj(Lj) - Bj-i(ij-i), (4.7) 

24 



where 



tti,2 = -iW ^(0,y)g J (W J (y),y;f,)dy (4-8) 

+Pj [ $ J s _ 1 (0,y)g J _ 1 (V^_ 1 (y),y;/i)dy (4.9) 

-Pj$?(0, L 3 )Pf (Lj ) Cj (Lj )ujj (4.10) 

'i'V (0, h)Pf (Lj)dj (Wj^^Lj) (4.11) 

-Pjfy-dO, -L j - 1 )P/_ 1 (^- 1 )c j _ 1 ( J L J -_ 1 )^_ 1 (4.12) 



-P i $3_ 1 (0,-L i _i)i^_ 1 (L j _ 1 )?l i _i(W i _ 1 ,W i _ 1 ;/i,L j _ 1 ). (4.13) 

By construction, the left hand side of (|4.7|) is invertible as a map from £^(0) x 

£"_ 1 (0) to £y. We will next estimate the terms in 72^2 and show that the right hand 
side is a perturbation. 

Note that for r = j — 1 , j the terms W r and W r depend on u> r as in Theorem 13.21 
In particular, the right hand side of (|4.7p depends on w|_x and Vj-i (through 
which involves as well as and u^+i (through which involves w") so that 

each equation is coupled to its left and right neighbour in the itinerary. 

Step 1: estimating P~5J) and flUD - 

We expand (|4.8[) and (|4.9[) . and determine the difference to AdjU- As a shorthand 
we use fj = f(qj\0). Upon expanding gj in w and /i the terms fj and d u fjWj cancel; 
similarly for gj-\- Thus, (|4.8[) and (|4.9[) . respectively, equal 

,0 

P 7 - / $«(0, y){d uu f 3 W] + d^f jf i)dy + 7^,3, (4.14) 

-Pi [ $U (o, iO(auu£#?-i + ^/»dy + ^- 4 , (4.i5) 

where 

^•3 = o(H 2 + Hll^-||+, i3 + (||Wi||+ iij ) 3 ) J 

K jA = G(H 2 + Wll^i-ill".^ + (||T^-_x||-_ 1> ^_ 1 ) 3 ). 
The estimates in Theorem 13.21 imply 

723,3 = 0(\n\ 2 + |^ s | 3 + R) + (|^ s | + P,)| M |), 
K jA = 0(\fi\ 2 + \Qj\ 3 + Pf_! + (|^| + Rj-iM). 

In the following, the estimates of Theorem 13.21 will be substituted directly without 
mentioning. 

We next estimate the term in (|4.14|) involving d uu fjWj. Since Wj is a fixed point 
of Qj^i we have 

Wj(y) = $5(|,, 0)wj + $;(y, L^PfiL^bjiLj) + H jtS ] + 72,- 6 . (4.16) 

Looking at the definition of Gj,i, the remainder term 72^6 consists of the integrals 
involving g 3 , while 72j,5 contains the terms involving Cj as well as £)j. 
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When estimating the order of TZja, (|2.2j) and the weight rjj allow to ignore the 
integrals and $j(x,Lj) in the sense that the integrals arc of order |<7j||^~. L .. This is 
estimated in (j5T81 to be of order \n\ + |u^| 2 + i? 2 . Concerning 1Z ji5 , (pH?)) and (jXTSj) 
show that TZ j>5 = 0{R 3 ). 

In order to get good estimates of the other terms in (|4.16[) we decompose 

duufj (Wj - '!-; (y, 0)^ + (y, 0)lu*) 2 = d /\ (W j - -!-; (y, 0)^ s ) 2 

+2a uu / J -[w j ,$5(j/,o)^] -^(^.(y.o)^) 2 , 

and substitute this into (jlTH)) . Consider X = J° $^(0, y)d uu f ] (W :j - $^(y, 0)wp 2 dy, 
and substitute (|4TTo) for W 3 -. Using (g^J), d3HJ) as well as TZjfi = 0(Rj) we find 

\X\ < C J\- K ^\d uu fj\ (c- K ^ L ^\Pf(L 3 )\R + \Kjfi\fdy 

r° 

< C( / e - " -'' dyR 2 + |7e,- 6 | 2 ) 

<C(i?f + |^| 2 + |^ s | 4 ). (4.18) 

The remaining terms from f|4. 1 7[) are directly estimated to be of the order \ujj\ 2 and 
WWjW+^l i.e., of order + + Rj). 

In summary, including 7^,3, the term (|4.8j) can be written as 
-P, f° ^(0,y)dJ jf idy + TZ^, 

where 

7^,7 = O (ImKImI + K s + Rj) + h s l 2 + Rj\Uj\ + Rj) ■ 

The completely analogous computation for Wj-i and (|4.9|) shows that (|4.9[) equals 

/ ^ 1 (0 1 »)8 |1 /j(id»t^, 

where 

Letting the bounds in these integrals tend to infinity generates an error of order 
(Rj + Rj-i)\fj,\. Hence, the sum of (|4.8j) and (|4.9p is of the form 

Mji ■ O R, - ■ 11,,). (4.19) 

Step 2: Estimate (|4~T0l) - (|4~T3l) . 

Using (|5T5|) and ([2~2]l shows that (|4TTT|) is of order 

R 2 j (\u]\ 2 + Rj\u s j\ 2 + RjR 2 + R 2 ) + Rj\fx\, (4.20) 
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and similarly (|3.16l) implies (|4.13[) is of order 



RU (|wS-i| 2 + ^-i + ^-i^-i + ^-i|^-i| 2 ) +Rj-iH (4.21) 

By ((33|) it follows that (|4TT0|) is of order 

i?f (l^l+^l), (4.22) 

and by (|3. 16[) it follows that f|4.12[) is of order 

^-iC^-il^il + lw^il). (4.23) 

In summary, lZj.2 = -Mjfi + 0(TZj^ + TZj^ + Rj\&j\ + Rj-i I)) anc ^ m particular, 
using (|4.4j) and (|4.5[) , the right hand side of (|4.7j) is of the form 

- Wj + O (N 2 + i? 2 + RU + | M | + K| 2 + |^_! | 2 )) (4.24) 

Step 3: Apply the implicit function theorem. 

On account of (|4.24[) the right hand side is a perturbation of the invertible left 
hand side for large Lj and and small \vj\. Notably, the constants in the above 

order estimates all depend only on C from Theorem 13. 21 and are uniform in Lj, 

Hence, for finite J and any closing condition there are E\,ei and L* such that 
the implicit function theorem immediately applies to the then finite system (|4.7[) for 
j G J° . This gives unique C k smooth solutions uij,u!^_ 1 and, due to (|4.24l) . there is 
C > such that these satisfy 

+ |^_! - Vj \ < C + \v 3 \ 2 + i? 2 + R^) . 

Since TZ is a function of u),v,L,fj, this proves items 1 and 2 of the lemma. Item 3 
follows from the above form of "R>j,2, and item 4 from this and the implicit function 
theorem. 

For infinite J the inverse of the map generating the left hand side of (|4.7p , 

n -> J] (^J + (wJ^i-i) ^ ^ - ^--i. 

is given componentwise by (Pj,(Pj — Id)) with norm measured by that of Pj := 
Proj(7J|, H^x). Since C ro( j is finite, this is uniformly bounded in j. The estimates 
of the right hand side of (|4.7p and that this involves only j and j — 1 immediately 
gives continuity in L,Q,v in the spaces £, f2, V. Smoothness of order k follows again 
using that (|4.7I) is local in the index j; smoothness in /z is straightforward. Hence, 
the implicit function theorem applies as in the finite case. □ 

4.2. Completing the solution using system parameters. Upon substitut- 
ing the solutions uip &J_i, j G J°, from Lemma T4.4I into the fixed points Wj,Wj of 
Theorem 13.21 the remaining parameters are v, fi and L. 

As explained in Remark HTTl the projection of (|4.1[) to the flow direction qj(0) 
is trivially solved. Therefore, the spaces Ej 1 C kerPj identify the remaining so- 
called bifurcation equations that determine /x via (|4.1[) multiplied by all a® r , j G J b , 
r G {1, . . . , dj}, which yields 

J2 <(W j -W f -i)(0;i;,L >A i),o5 i X r = 0. (4.25) 

r=l,...,dj 
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Recall that a° r form a basis of so that each summand has to vanish. 

We will show that the boundary terms bj enter the j-th bifurcation equation via 

4>(L) := (Bj(L) , a° +l r ), 
Bj,r{L) := {BJ(L) , a°, r ), 

Bj(Lj-l,Lj) = J2 (4-l,r(ii-l)+B J >(i J 0)o?,r- 
r=l,...,dj 

To capture the leading order dependence of the j-th bifurcation equation on the 
neighbouring (j ± l)-st we define (recall P, from the proof of Lemma EQ)) 

G, := P i+ ife+i + B» +1 (L J+1 ) + Bj(Lj)) 

Gj-l := Vj- X + (Id - /' + ^(L,-.!)), 

S+. := -(-i>](0,L J )P^L J )^(-L J ,0)G J ,al r )dy, 

Sj,r-= ($;- 1 (0,-^-i)F;_ 1 (L J -i)^_ 1 (i,_ 1 ,0)G J _ 1 ,aO r )dy, 

5 i = Sj(vj+i,v j -i,L j ,L j -- i ,L j+1 ) := ^ (5^ r + 5+ r )a° r . 

r=l,...,dj 

Note how in these terms the coordinate parameters Vj±i arc transported by the linear 
evolution from <?j±i(0) to ?j(0), respectively. 
The estimates ([2"^2"jl . (|3~3]) and (gTSJ imply 

^'(Aj-ii-kj) = 0(Rj_i + Rj), (4.26) 
St - 1 + i?, 2 +1 + (4.27) 

Sr. = O^d^-il + Pf + Pf.J). (4.28) 

In the following lemma a subtlety for the correction Sj in the error terms Rj, Rj 
arises. For simple leading eigenvalues Sj = is possible everywhere, except in the 
estimates of Theorem 13.21 for j € Jp, which requires an exponential weight r/j > 0. So 
far this was irrelevant, but now it becomes important, and therefore we let Rj = Rj, 

Rj = Rj denote error terms where Sj can be set to zero for simple leading eigenvalues. 
Lemma 4.5. 

1. For j G J b and L, v as in Lemma \4-4[ equation (|4.25[) is of the form 

■MjV = Tjivj) - Bj(Lj-i,Lj) + Sj(vj+i,Vj-i, Lj, Lj-!, Lj+i) + Rjs + Rj,w 
where the remainder terms are C k smooth and 
Kj, 9 = o(\ii\(Rj+Rj- 1 + \v j \ + \ti\ 



Kj,w = o (M (N + RjRj + Rj-iRi-i 

+ \vi+i\R 2 j{Rj + Rj+i + \v j+ i\) + {vj-xl^^Rj-i + Rj-i + \vj-tl) 
+Rj(Rj + Rj + Rj+i) +Rj-i(Rj-i + Rj-i + Rj-2) 

2. Under Hypothesis^ for sufficiently small \v\ and large L, there exist 
unique C k smooth Rj,n of the same order as Rj,io such that (|4.25|) is solved 
if and only if for j G J h the itinerary parameters fj,* satisfy 

fij = Tj(vj) - Bj{Lj- U Lj)+Sj{vj +1 ,Vj-i,Lj,Lj_ 1 ,Lj +1 ) + Kj,n- 
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Here, M.j ,Tj,Bj, Sj have reducible indexing (they do not differ on repeated parts 
of the heteroclinic network in the itinerary) . 

Proof. Substituting the definition of Qj, and suppressing some variables for read- 
ability, the summands in (|4.25[) are 

- ($?(0, L^PJiL^ + $^(0, -L J -_ 1 )P?_ 1 (L,-_ 1 )6 i _i,a° r ) (4.29) 

-<#j?(0 > L i )i^ , (i J -)c i ,o5 ir > (4.30) 

-(<S>](0,L 3 )P»(L^,al r ) (4.31) 

-($5_ 1 (0,-i 3 -_ 1 )J^_ 1 (L i _i)c,-_i,o5 ir > (4.32) 

-<$5_ 1 (0 > -L J _ 1 )J^_ 1 (i i _ 1 )d 3 -_ 1 ,o5 >r > (4.33) 
,o 

+ (/ ^(0,y) ft (W 3 ,2/;/i)dy,a° r } (4.34) 

-(/ *5_i(0,»)S J _i(W i _i,y;^)dy,o5 ir > (4.35) 

Note first that term (|4.29p is precisely — B 3 -_i r (ij-_i) — Bj iT (Lj). 

Step 1: Estimate l|£30 )H|O5p . 

From the estimate (|4.20[) of (|4.11[) and Lemma |4~47 2) we infer that (|4.31|) is of 
the order 

R) (k) + |u 3+1 | 2 + R) +1 + + 4-(H 2 + i?| + fl|_ x )) + J? 3 |/i|. 
Similarly, now using (|4.21[) , (|4.33|) is of theorder 

R 2 j-i (^ 2 -i + l^-il 2 + ^ 4 -2 + Rj-i&j-i + Rj-i(\v 3 \ 2 + + RU)) + 

Substituting the expansion of from Lemma l4"4T 3) into (|4.30p and using (|3.9j) 
as well as Lemma [474^2) for \uij\ gives 

s+ r + o (fi (| M | + |^ + i,i I + Rj(\vj\ + R) + ^ 2 -i))) • 

Similarly, using the expansion of £^_ l5 (14. 32[) contains S~ r , and, by (|3.10[) . the rest is 
of the order 

4 2 _ x (| M | + |fti- M | + Rj-iQvjl + R) + RU)) . 

In summary, after some computation, (|4.30p + (f4.3ip + (|4.32| + ()4.33p minus S- r is of 
the order 

\ti\(Rj + Rj-i) + IvjKRtRj + R^Rj-i) 
+\v j+1 \R 2 j (R j + R 3+1 + \v j+1 \) + Iv^lR^iRj-i + Rj-i + I) (4.36) 

■R?i' : ; + R) + R"+i) + 4 2 -i(4-i + R "-i + R *-2)- 

Note how remainder terms come from the local piece of the itinerary, but also from 
one and two steps further along the itinerary, if the itinerary is that long. If the 
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itinerary is shorter and not 'per' then all the terms with indices outside the range of 
the itinerary vanish by definition. 

Step 2: Expand and estimate (|4.34j) and (|4.35|) . 

Similar to the proof of Lemma l4~4l the idea is to expand (|4.34[) and (|4.35[) . so that 
the sum can be written as 



-Mjn + Tjiv^+Kj,!!. 
We write fj = f(qf, 0) and expand g 3 - and g -\ so that (|4.34j) + (|4.35|) equals 

< / $"(0, y)(d uu fjW? + drfjridy, oL) (4.37) 



j 

Li 



where 



$^(0, yKduufjW?., + dufjridy, a° r ) + K j>12 , (4.38) 

-Lj—i 



+(II^IU j ,l j ) 3 + (||^iIU j _ 1 ^_ 1 ) 3 . 



Theorem 13.21 and Lemma 14.41 imply 

||^(a,>,^)lk^ = OQvl\ + Ivjl+Rj + R 2 ^), 
WWi-iiuf^^Lj-^W^L^ = 0(\fx\ + \v j \ + R 2 + R j - 1 ), 

which yields 

Hj,u = 0(ImI 2 + N 3 + R* + Rj-i + (Rj + Rj-i + \v 3 \M). 

In (|4.37[) we write 

duufjWf = duufjlWj - $5(x,0)w3] 2 + d uu f 3 (^(x,0)iu s 3 ) 2 



(4.39) 



and consider the resulting integrals in (|4.37p from each of the three summands (I-HI) 
of the right hand side of (|4.39[) . Since we already have an error of order of ^,12, we 
focus on the additional contributions. 

(I) The integral over the first summand gives rise to the term X from the proof of 
Lemma l4~4l Using (|4. 18[) it does not contribute further to the order of Hj,i 2 - 

(II) The integral over the second summand d uu fj(^ S j(x, 0)wj) 2 reads 

(d uu f j ($ s j (y,0)v j ) 2 ,a j , r (y))dy + 0(|££|(KI + |^ s |)). 



Substituting the estimate of Lemma 14.41 implies that the remainder term of this is of 
order 

(|/i| + \ Vj \ 2 + R 2 + R 2 ^)\v 3 \ + \fi\ 2 + R$ + < C\K jA2 \. 
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(Ill) Using (|4.16|) and the same integral computation as for the estimate of X in 
Lemma \A. 41 the integral over the third summand d uu fj[Wj — ^(x,0)uij,^(x,0)u)j] 
is of order iT^j^Hw^l with TZj t e from that proof. Lemma 3^2) shows 

Hjfi = O (| M | + K'jf + + Rf) = 0(TZ jll2 ). 

The same holds for the corresponding terms in (|4.35[) . Hence, going to infinite 
integral bounds as in the proof of Lemma l4~4l we obtain that (|4.34[) + (|4.35|) can be 
written as 

Combining this with the remainder term (|4.36[) of step 1 proves part (1) of the lemma 
statement. 

Step 3: Apply the implicit function theorem. 

Recall the change of parameters in the Lyapunov-Schmidt reduction of Aifi = X 
discussed after the definition of M in (|4.3[) . The itinerary parameters are such that 
Mjlpj : M. dj — > Ej is the identity (on the chosen basis) and generates the solvability 
conditions Xj = Xy . Using part 1 of the lemma, estimates (|4.26p , (|4.27|) and (|4.28|) 
and Hypothesis [4] allow to apply the implicit function theorem immediately for finite 
J. 

As in Lemma 14.41 all estimates are local in the itinerary, i.e., in the index j, so 
that smoothness in (v,L) £ V x C(L*) follow. Hence, the implicit function theorem 
also applies for infinite </, which proves part (2) of the lemma statement. □ 

The basis to study the leading order geometry of bifurcation sets is the following. 

Lemma 4.6. There exists L* > L 1 such that the following holds. 

1. Assume Hypothesis^ with i>j = — K s + icr|. There exist /3j^ r ,hj. r £ R, /3| £ R™, 
hj £ Ej such that for all L £ Kj(L*) and <5^ < min{Kj, p)} 

B S 3 {L) = e / C,»(2^/. + fifth) + C(c" (2K H' 5 P L ), 
B jtr (L) = e~ 2 ^ L cos(2ct j s L + fc, r )hj, r + C( C - (2 ^ + ^ )L ). 

If Hypothesis \3{2) and\9H) hold, then there is r s £ {1, ...,dj} such that 
h j>r s + 0. If Hypotheses^ 2) and\9(6) hold, then h] ^ 0. 

2. Assume Hypothesis^ with vj = -K u ±ier". There exist $^ r , hj. r £ R, /3" £ R n , 
h s j £ Ej such that for all L £ Kj(L*) and <5" < min{Kj,p"} 

BjiL) = c' 2 ^ L Cos(2a j L + /3%)h] + 0( C -^ +s ^ L ), 
B jtr (L) = e~ 2K 1 L cos(2ajL + $^ r )h^ r + 0(e-^ +s ^ L ). 

If Hypotheses Wfl) and\^3) hold (replacing j — 1 by j), then there is r u £ 
{1, ...,dj} such that hj—i r u ^ 0. If Hypothesis and (5) hold, then 

i>'; ■■ «■ 

Proof. The expansions are a direct consequence of substituting Lcmma [2.3f 4) into 
the definitions of B j<r (L), B^ r (L), B S / U (L), using Lemma [2~3l 2). and angle addition 
formulae. □ 

Proof of Theorem 14.31 Concerning the expansion of fi*, we substitute the 
results of Lemma 14.61 into the expansion in Lemma I4.5f 2) and solve the resulting 
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equation for \i. The exponents Sj in the remainder term are not restricted further 
than in Lemma 14.61 since the relevant terms in TZj^o of Lemma 14.51 are i?| and Rj-\- 
The terms from are, using placeholder terms Y, Y, of the form 

{^L^PfiL^i-L^^Y^lr) 

($»_ 1 (o,-i j _ 1 )if_ 1 (L i - 1 )$»_ 1 (i j _ 1 ,o)y, a ° r ). 

As noted in the proof of Lemma |2~37 4) we can replace ^_ 1 (x, 0), <&^(x, 0), 0), 
$"(a;, 0) by the evolution and projections of the variation about Pj to leading order. 
These linear evolutions generate terms of the same form as the leading order part of 
the expansion of the nonlinear terms in Lemma 14.61 

The cosine terms stem from resolving Fj, Fj-\ in terms of the leading eigen- 
values, which gives sine and cosine terms with the same angle arguments appear 
with Vj±i -dependent coefficients, respectively. After multiplication with MJ 1 in each 
component all terms can be joined into one cosine term for and one for with 
Wj-ti-dependent phase shift, respectively. This gives the diagonal nature of the linear 
coefficient maps. Expanding phase shift and cosine, only the linear Wj±i-dependence 
is kept and the rest subsumed into the remainder terms. 

Finally, the order of -^-Tij follows from inspecting (|4.30|) - (|4.35|) as in the proof 
of Lemma [4.51 This completes the proof of Part 1, and Part 2 follows from this in 
combination with the estimates of Theorem 13.21 and Lemmata 14.41 14.51 

Part 3 of this theorem is a consequence of the fact that fixed points of Qj arc coor- 
dinates of trajectories, see Corollarv l3.3t and the uniqueness of solutions in Lcmma [4.4l 
and 14.51 Injcctivity from (u, L) to trajectories for fixed fi, follows the expansion in 
Lemma I4.4f 3) and Corollary I3.3( injcctivity in p follows from invcrtibility of A4 ro d- 
It remains to argue that this covers all solutions in a neighbourhood of the visited 
heteroclinic connections. 

The assumptions on Xj guarantee that the j-variations obtained from u are co- 
herent with C and so small that Theorem 13.21 and Lemmas 14. 4[ 14.51 apply. Since fixed 
points of Qj are surjective onto such j-variations in a neighbourhood of each hete- 
roclinic connection, each of these corresponds to a unique fixed point of Qj near pj. 
Hence, the coordinate parameters must solve (|4.1[) for the given /i. The Lyapunov- 
Schmidt reduction of this system done in Lemmas 14.41 14.51 provides necessary and 
sufficient conditions on all such solutions. Therefore, the coordinates parameters v 
and travel times L derived from u must solve part (1) of Theorem 14.31 In particular, 
these solutions indeed cover all trajectories of that remain in a neighbourhood 
of the selected itinerary. 

5. Sample bifurcation analyses. If the itinerary C does not have repeated 
elements of C^, then p, = /i* so there arc no solvability conditions for (v,L) in Theo- 
rcm l4.3l Hence, it already proves existence, uniqueness, and the parameter expansions 
of certain periodic, homoclinic and heteroclinic orbits as follows. If the heteroclinic 
network allows for recurrence, e.g., if p\ = p m , these solutions are simply recurrent, 
or 1-recurrcnt, in the sense that during their minimal (possibly infinite) period they 
intersect each section £,,- at <?j(0) for all qj £ C2 precisely once. We call such solutions 
1-periodic, 1-homoclinic or 1-heteroclinic orbits. 

Corollary 5.1. Assume Hypotheses [7] anrf[^} // J° = J rc d (in particular J = 
{1, . . . ,m} is finite and J ro d unique) then there exists s > and a neighbourhood IA 
of the itinerary such that the following holds. 
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1. If pi 7^ p m then the set of 1-heteroclinic orbits from pi to p m in IA is bijec- 
tive to the non-empty set of fixed points of Lemma \4-4\ an d parameters as in 
Lemma \4-5\ under the 'het' closing conditions L\ = L m = oo. 

2. If pi = p rn then the same holds for the set of 1-homoclinic and 1-periodic 
orbits under the 'horn' (L\ = L rn = oo) or 'per' (L\ = L m ) closing conditions, 
respectively. 

3. Under the assumptions of Theorem \4-3^ 1) in each case v and L parametrise 
the corresponding solution set, and p, G M. d satisfy the given expansion. 

In the following we illustrate for some basic heteroclinic networks how to deter- 
mine and use the equations for of Theorem l4.3l Most of the results are well known, 
but we hope to show the ease in obtaining them using that result. The last class of 
examples are for heteroclinic cycles between one equilibrium and one periodic orbit 
mentioned in the introduction. 

5.1. Heteroclinic orbits. The simplest heteroclinic network consists of two el- 
ements pi ^ p 2 connected by one heteroclinic orbit q 2 , which enforces the 'het' closing 
condition. We thus investigate the existence and variation of the heteroclinic connec- 
tion q2 upon parameter changes. Note that L\ = L2 = 00 so that the bifurcation 
equations only contain terms of order \/i2\ e , \v2\ e+1 , £ > 1. 

In this context E 2 = H f + C £1(0) + -^1 (0)> where the inclusion is due to the 
flow direction, which was removed for equilibria. The number of bifurcation equations 
is c^2 = dim £2 = n — dim(iJf + H\ l ) — 1, and the number of remaining coordinate 
parameters V2 for potential tangent directions is given by dim(_ff| n Hi). 

Typically, dim(i52) is maximal, so that no tangencies occur and is minimal, 
which implies = n + 1 — dim(W u (pi)) — dim(W s (p2)) is the codimension of the 
heteroclinic connection and #J* = dim(W u (pi) fl W s (p2)) — 1 is the dimension of the 
set of heteroclinic trajectories. 

5.1.1. Saddle-saddle connection. Suppose both p\ and P2 are equilibria or 
periodic orbits connected in a saddle-saddle situation dim(W u (pi))+dim(W s (p2)) = n 
for which the linear codimension is 1, i.e., > 1 and typically d% = 1. Since 
the heteroclinic connection cannot be transverse dimRgPj counts tangent directions 
(except the flow) in W u (pi) D W s (p2)) and typically is zero. 

In the typical case the bifurcation equation from Lemma [4.5l reads M.2^2 = C(m!)> 
[12 G K, and if M.2 ^ a heteroclinic connection exists only at /i2 = 0. If there is an 
auxiliary parameter fi2 we obtain a second contribution to the Mclnikov integral and 
the leading order bifurcation equation 

M2H2 + Mifi-2 = 0. 

Hence, if A^2 and M2 are non-zero, then heteroclinic orbits exist locally on a curve 
in the parameter plane. 

5.1.2. Tangent source-sink connection. In a source-sink case, where the het- 
eroclinic connection has non-positive linear codimension, generically = 0, i.e., there 
are no bifurcation equations for parameters. In that case, the coordinate parameters 
for heteroclinic orbits are given in Lemma 14.41 Note, that for negative linear codi- 
mension coordinate parameters Vj appear. 

If the heteroclinic connection is tangent, i.e., d2 > 1 we have J' = {2} and for 
single tangent direction 4 = 1| "2,^2 G R- The quadratic function 72(^2) is then 
scalar and can be written as av\ for a G K so that the bifurcation equation reads 

M2V2 =avl + 0(\v2\ 3 + (4)- 
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Hence, heteroclinic connections typically (a 7^ 0) occur on a parabola in the (^2,^2)- 
parameter plane at leading order. 

Including an auxiliary parameter JI2 we can trace tangent heteroclinic connections 
in the (/X2, /^-parameter plane. Tangencies are located at roots of the derivative of 
the bifurcation equation with respect to w 2 ■ At leading order this gives V2 = so that 
to leading order tangent heteroclinic orbits lie at M2H2 + M-2p-2 = 0. 

5.2. Bifurcations from homoclinic orbits. The situation of the previous sec- 
tion for pi = P2 allows for more interesting solutions and the 'horn' as well as 'per' 
closing conditions. To serve readability we omit the subscript which labels the single 
equilibrium or periodic orbits in the following. 

We consider the generic transverse case H s n H u = {0} where no parameter i>2 
occurs. Since d,2 — 1 in that case, all reduced index sets contain only one element, 
and for those we omit the labels. Hence, for any itinerary, each of the bifurcation 
equations reads, with fx, ft, 7 G R, 

fi = c ~ 2kUl j cos(2a u Lj + /3)C + c" 2 ^^- 1 cos(2cr s L J _ 1 + 7 )£. (5.1) 

Here we omitted the term IZj and set = £'j = since these terms do not appear at 
leading order in the following considerations. The occurrence of parameters L r and 
the range of indices depends on the choice of itinerary. 

We first consider an equilibrium p where C(L*) is continuous, and then a periodic 
orbit p where C(L*) is discrete (and the un/stable dimensions change). 

5.2.1. Equilibrium at p. The analysis of homoclinic orbits that do not pass by 
the equilibrium p, i.e., of 1-homoclinic orbits, is the same as in Section [5.1. H 

2-homoclinic orbits. Homoclinic orbits could pass by p any number of times before 
connecting to the un/stable manifold. Striving for illustration, we only consider 2- 
homoclinic orbits that pass by p once. 

The itinerary C for these orbits is as in Figure ITTT1 in |JT]and contains three equilib- 
ria pj = Pi, j = 1, 2, 3, under the 'horn' closing condition. The un/stable dimensions 
are all the same, respectively, so d,2 = dz = 1. When applying Theorem 14.31 with 
the 'horn' closing condition, one free parameter L2 appears (note L\ = L3 = 00), see 
Figure ll.l| and from (|5.1[) we obtain the system of bifurcation equations 

H = c~ 2kUL C cos(2cr u i + 0) 
M = e - 2 « Si £cos(2cr 8 £ + 7), 

the first for j = 2 and the second for j = 3. Equating the right hand sides gives the 
solvability condition 

c- 2KUL Ccos(2cr u L + P)= e- 2KSL £cos(2a s L + 7). 

If k s k u , say k u < k s , the leading order equation is 

e- 2K " L C cos(2ct u L + 13) =0, 

Hypothesis [5] implies ( ^ so that solutions at leading order exist if and only if 
<T U 7^ 0, which is the well-known Shil'nikov saddle-focus configuration. The arising 
infinite sequence of solutions persists (due to transvcrsality) under the higher order 
perturbation of the remainder term. 
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Concerning vanishing coefficients of leading terms, we outline the result mentioned 
in Remark 11.11 in case 1Z is higher order with respect to the terms in (|5.1[) (roughly 
speaking this is valid for small difference of leading stable and unstable rate, and large 
gap to the next leading rates). We thus keep all terms from (|5.ip and obtain 

e 2L (K "- K =) = Ccos(2q"L + /3) 

£ cos(2cr s L + 7) 1 ' ' 

For real leading eigenvalues a sign change of Q for £ 7^ implies the emergence 
of a solution from L = 00, i.e., for sufficiently large L for Theorem 14.31 to apply. 
As mentioned, the results in [33] provide a complete study of this situation and, in 
particular, do not require such restrictive spectral configurations. 

In the resonant situation k s = k u equation (|5.2[) applies as well with left hand 
side equal to 1. This yields solutions if either a s = or a u = 0, or else under non- 
resonance conditions on these and 7, j3. Typically, infinitely many solutions persist 
when including the higher order terms. For real leading eigenvalues there is no solution 

1- and 2-periodic orbits. Periodic orbits could pass by p any number of times 
each period, but here we only consider the cases where this number is one or two. 
Typically stable and unstable rates differ, say k u < k s , and the coefficient £ 7^ 0. For 
the 1-periodic orbits the leading order equation according to (|5.1[) is 

fj, = c- 2 ^ L £cos(2cr s J L + 7). 

The well-known result follows that only for cr s 7^ periodic orbits at fj, = co-exist 
with the homoclinic orbit, and in fact accumulate on it. 

The 2-periodic orbits are encoded in the itinerary of the 2-homoclinic orbits with 
'per' closing conditions so that indices of L need to be taken mod 2 + 1 and (|5.1|) 
yields 

/1 = c- 2KUL2 Ccos(2cr u L 2 + /3)+ e- 2KSil £cos(2cr s Li + 7) 
H = c-- 2kULi Ccos(2ct u Li + 0) + e- 2KSi2 £cos(2cr s L 2 + 7), 

with the two parameters L2 and L±. Both equations coincide with that for 1-periodic 
orbits if L2 = L\. 

In case k s > k u the solvability condition to leading order as L* —J- 00 is 
e - 2KUi2 Ccos(2cr u L 2 +/?) = e- 2KUil Ccos(2cr u L 1 + /?), 
or equivalently, for ( 0, 

^( Ll -L 2 ) _ COS^"^ + 0) 



c 



cos(2cr u L 2 + 0) ' 



which shows that for the Shil'nikov saddle-focus with cr u ^ there are infinitely many 
persistent solutions. 

Assuming that It is higher order with respect to all terms in (|5.1[) (i.e., the 
spectral configuration is as mentioned in the 2-homoclinic case) the resulting leading 
order solvability condition for cr u = <r s = and £ 7^ 0, can be written as 

-2k s L 2 _ p -2k s L! f 

(5.3) 



-2k u L 2 _ c -2k u Li £ 
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We observe that a sign change of £ leads to the period doubling bifurcation of a 
solution curve L 2 <~ L\ with either L 2 > L\ or vice versa. Note the analogy to the 
2-homoclinic case. 

In the resonant case k s = k u , for cr s / u = the left hand side of (|5.3[) is 1 so that 
solutions do not exist for £ ^ £. 

5.2.2. Periodic orbit at p. The form of the abstract bifurcation equations does 
not change much when p is a periodic orbit, only L is discrete, but the set of solutions 
near the homoclinic orbit to p may change dramatically. 

The reason is that such a homoclinic orbit is generically codimension-0 since 
the flow direction is the center direction which counts towards stable and unstable 
dimensions. Indeed, typically the complement to H u + H s only contains the flow 
direction so that d = 0. Hence, there is no parameter needed and no solvability 
condition. This means that under the assumptions of Lemma 14.41 solutions for all 
itineraries and for any small parameter perturbation can be constructed by adjusting 
the coordinate parameters (w s ,w u ) according to the expansion in that lemma; note 
that all constants depend only on p and q, and in particular are uniform for all C. 

Complicated dynamics typically occurs since the diffeomorphism generated by any 
suitable Poincare map has a transverse homoclinic orbit which is one of the paradigms 
of chaotic dynamics |30j . Note that in the present setup the ambient dimension is 
arbitrary. 

We next show how Theorem l4.3l can be used to prove conjugacy of the dynamics 
on the local invariant set to shift dynamics on two symbols. Let y C U be the set 
of trajectories {u(x) : x £ R} contained entirely in the neighbourhood IA of q from 
Theorem I4.3T 3). Take a suitably large Poincare section E p transverse to the flow 
containing p(0). For all solutions in y with u(0) £ Si (without loss of generality), we 
find a unique sequence x s £ R, s £ Z such that xq = 0, x s < x s+ i and x s+ \ — x s is 
minimal so that u(x s ) 6 Ei U S p for all s. 

This defines a unique symbol sequence (a s ) se z, a s £ {X, Y} by setting a s = Y 
if u(x s ) £ Si, and a s = X if u(x s ) £ E p . Since we require a minimum travel time 
L* from E to E (the time from E to E is constant) there is a well defined minimal 
number jy{L*) of X's after each Y in the sequence. 

Corollary 5.2. Assume (jl.ip possesses a homoclinic orbit q to a hyperbolic 
periodic orbit p, and suppose that dim(if u + H s ) = n — 1. Then there is a number 
L* > and a neighbourhood U of q such that the invariant set y in hi is bijective 
to the set of sequences {(a s ) se z} defined above for which there are at least jy(L*) 
symbols X after each Y . 

Proof. In this case there is no system parameter \x and no coordinate parameter v 
appears in Theorem l4.3l Hence, there is a minimal travel time L* and a neighbourhood 
Li of q such that the solutions for all itineraries in that neighbourhood are bijective to 
the sequences of travel times in C{L*). Since any orbit that lies in U for all time has 
a unique such sequence, the entire invariant set y in U is bijective to the sequences in 

C* := {C(L*) : C is an itinerary}. 

In particular, any orbit in y has a unique itinerary of intersections with E p and Ei as 
defined above, i.e., the map from travel time to these symbol sequences is injective. 

To prove surjectivity consider a symbol sequence {(a s ) sg z} with at least jy{L*) 
symbols X after each Y . We construct an itinerary that generates a solution with 
that sequence. For this we define a subsequence (& s ) se s, B C Z of a s and then 
consider the itinerary generated by the sequence of Y's in b s . First remove jy(L*) 
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symbols X after each zero in (a s ). If the resulting sequence is periodic, then b s is 
a minimally periodic subsequence, say of period m, and we employ the 'per' closing 
conditions. If the resulting sequence is constant for s > s + and/or for s < s_ then 
(these must be constant X's) b s is defined as the sequence between such maximally 
chosen s_ and/or minimally chosen s + . If the resulting B is bounded we employ the 
'horn' closing conditions, else the corresponding 'semi' condition. 

Let C be the itinerary pj = and qj = q*, where J is bijective to {s : b s = X} 
and min J = 1 if it exists. There is a unique solution obtained by Theorem 14.31 for 
that itinerary with Lj = IfT where lj > jy{L*) is the number of consecutive X's in 
b s that follow the Y corresponding to j in the bijection between J and {s : b s = X}. 
The above defined map from y to these symbol sequences thus surjective. □ 

Corollary 5.3. The dynamics of on the invariant set near a trans- 

verse homoclinic orbit to a hyperbolic periodic orbit (i.e., the trichotomy spaces satisfy 
dim(i7 u + H s ) = n — I) is conjugate to (suspended) shift dynamics on two symbols. 

Proof. Let L* be as in Theorem 14.31 in this setting and y the local invariant set. 
Let (c r ) rS z be a bi-infinitc sequence of symbols and 1. We define the sequence a s 
with symbols and meaning as above: replace all 1 by Y followed by jy{L*) X's, and 
replace all by X. 

By Corollary 15.21 there exists a unique orbit solving (jl.lj) corresponding to that 
sequence. By construction, there is a unique sequence of time steps x r , r £ Z such 
that the time evolution of the trajectory u for the unique discrete times of intersection 
with Si and E p is precisely the shift of the index c r i— > c r +\. 

Concerning continuity of this map from trajectories to symbol sequences, we 
consider the usual product topology on symbol sequences where cylinders are open 
sets, i.e., sets with some prescribed finite sequence of adjacent symbols. The norm 
generating this topology is given by (|4.6|) when taking Lj £ {0, 1}, see, e.g., [To] . 

In the direction from y to symbol sequences continuity follows from the con- 
struction: convergence in y means that initial conditions converge, which implies 
convergence of travel time sequences L in C(L*). By construction of the symbol 
sequences this implies their convergence in cylinders. 

Conversely, let a and a' be symbol sequences so that a — > a' in the cylinder topol- 
ogy. By construction, the travel time sequences L and L' of the corresponding solu- 
tions in y converge in £(£*). Since the coordinate parameters Qj and are con- 
tinuous in L this implies that the coordinate parameters converge as in Lemma l4.4f 2). 
Therefore, the solutions in y converge as well. □ 

If the homoclinic orbit q\ is tangent, e.g., dim(iT u + H s ) = n — 2, the bifurcation 
equation for 1-homoclinic orbits is the same as for the tangent source-sink hcteroclinic 
in Section [5.1.2l The dynamics near such a homoclinic tangency is more complicated 
than in the above case, see, e.g., [30] . 

5.3. Bifurcations from EP heteroclinic cycles. In this final section we con- 
sider heteroclinic cycles between one equilibrium p\ = E and one periodic orbit P2 = P 
with heteroclinic connections qEP = qi from E to P and qpE = qi from P to E. Such 
cycles have been recently found in a number of models, see [H 21 |3H I3H] and the 
references therein. EP cycles are also called singular cycles, and have been studied 
from an ergodic theory point of view in [TJ [23l [24] [27J [29] , and further papers by these 
authors, looking for instance at properties of non- wandering sets. 

Generally, in an EP cycle one connection is generically transverse, while the other 
has a positive codimension, see [32] . Here we consider the following three cases: 
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Fig. 5.1. Schematic plot of a segment of a general itinerary in an EP1 or EP2 heteroclinic 
cycle. Passing through T,ep does not require a parameter fJ.j±i since the heteroclinic connection 
from E to P is transverse, but in the EP2 case a parameter Vj±\ S R appears. On the other hand, 
passing through Epg at position j generates itinerary parameters fij £ R d for the EPd case, since 
this connection has codimension d = 1 or d = 2. 



EP1: the connection from E to P is transverse and one-dimensional, and the 

connection from P to E is codimcnsion-1, 
EP2: the connection from E to P is transverse and two-dimensional and the 

connections from P to E is codimension-2, 
EPlt: the connection from E to P is codimcnsion-1 and the connection from P 

to E is tangent. 

Concerning stable and unstable dimensions at E and P let ie be the number of 
unstable dimensions at E and ip the number of unstable dimensions at P including 
the flow direction. Let <1ep = di be the codimension of the heteroclinic connection 
from E to P and dpE = d\ that for the connection from P to E. The three cases are 
then as follows: 

EP1: i E = i Pl d PE = 1, d EP = 0, 

EP2: ip=i E -l,dpE = 2,d E p = -l, 

EPlt: ie = ip — 1, dps = d E p = 1- 

As in (|5.1|) for the homoclinic orbits to an equilibrium, we use this and Theo- 
rem 14.31 to obtain the form of the bifurcation equations without choosing a specific 
itinerary now, and omitting IZj and BJ +1 , B^_ 2 for the moment. The occurrence of 
parameters Vj , Lj and the range of indices then depends on the choice of itinerary. 

For the EP1 case only the connection from P to E generates a parameter, i.e., if 
index j corresponds to this connection then there is no bifurcation equation for indices 
j ± 1 for any itinerary, see Figure 15.11 Therefore, we denote the leading eigenvalues 
with subindices E and P as they always appear in the same order in the bifurcation 
equations which have the form 

H = c - 2k z l ' cos(2a E L j + /3)C + e" 2 ^^'- 1 cos(2oj,i J -_i + (5.4) 

For the EP2 case the connection from E to P is transverse and hence does not 
require system parameters. But the set of heteroclinic points from E to P is two- 
dimensional, i.e., a parameter Vf € R arises for j' corresponding to that connection. 
In any itinerary the index j = j' ± 1 then corresponds to the connection from P to 
E, which has codimension-2 so that /i* from Theorem 14.31 is two-dimensional and we 
write fXj = (//i,/^) since this itinerary parameter always corresponds to the same 
system parameters. See Figure 15.11 for illustration. Note that in this case £, £, and 
the image of are diagonal 2-by-2 matrices. With subindices E and P as for the 
EP1 case, the bifurcation equations read, for r = 1, 2, 

Hr = c- 2k e l i cos{2a E L j + f3 r + (3' r v j+1 )(C r + Oj+i) 

+e -2«5,L 3 _ 1 cos(2c7 u L . +lr + 7 > J _ 1 )(e r + (5-5) 
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(a) (b) 

Fig. 5.2. Sketch of an itinerary for a 2-homoclinic orbit for an EP1 heteroclinic cycle, (a) the 
heteroclinic cycle composed of q^, q| (solid), and a 2-homoclinic orbit (dashed), (b) Schematic plot 
of the itinerary (solid) for a 2-homoclinic orbit (dashed) for which L\ = L5 = 00. Note that the 
passing through £1 does not generate a parameter since the heteroclinic connection from E to P is 
transverse. However, each passing through £2 generates a single itinerary parameter. 



where Lj G [L* , 00) measures the time spent between Sp^ and £pp, while the discrete 
Lj-i G Kp(L*) approximately measures that between E^p and £p_E, see Figure [STTI 

The parameters Vj±\ € R can be viewed as varying the underlying reference 
heteroclinic connection from E to P. If the heteroclinic set from E to P has a winding 
number this can be used to obtain a continuous parameter Lj—i for the travel time 
near P, see [32] . 

In the EPlt case all /i* G R are one-dimensional and the connection from E to P is 
codimcnsion 1 while the connection from P to E is tangent. The tangency generates 
coordinate parameters Vj±i G R and we write 72(f) = av 2 for a G R. Note that 
Vj±i also appear in the bifurcation equation for the codimension-1 connection which 
is neighbouring this in any itinerary. Due to this coupling the resulting bifurcation 
equations cannot be reduced to the same basic building block form above, but to 

M2 = av)_ x + c^bL.-i C os(2(t|L,_ 1 + /3 2 )C 2 + «' 2 " ]J ' ' cos(2cr p i J _ 2 + 72 )£ 2 
,ui = c - 2k p l J cos{2a u P L ] + + /3[v j+1 )(Ci + ([v j+1 ) 

+e -2«» E L 3 -_ 1 cos ( 2o ^i J ._ 1 + 7l + 7 ^ J -_ 1 )(Ci + &V3-1), (5.6) 
M2 = av 2 +1 + c~ 2k bL ]+ i cos(2a^L i+ i + ^ 2 )C2 + c- 2K r L i cos(2a P L :j + 72 )6- 

Here /ii unfolds the connection from E to P and /i 2 from P to P. The first and last 
equation form a solvability condition, if the itinerary under consideration is that long. 

5.3.1. EP1 and EP2. The loci of simply recurrent solutions whose itinerary has 
no repetitions are explicitly given for the EP1 and EP2 case by the above equations: 
for 'horn' set either Lj = 00 or P/-i = 00 and note that Lj-\ is discrete; for 'per' any 
(Lj, Lj—i) for j = 1 with Lj_\ discrete generates a solution. 

For illustration we next consider 2-homoclinic orbits in the EP1 case; this also 
indicates complications arising from the terms omitted in (15.41) . 

2-homoclinic orbits to E pass by E once and P twice so that the itinerary is as 
in Figure HT21 and gives three parameters L 2 , L^,L^ where L 2 ,L4 are discrete. Since 
in this case we will get a solvability condition where the error terms have exponents 
from different pj, they need to be treated more carefully. From (|5.4[) the bifurcation 
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equations including IZj and £? E , Bp as error terms, are 



H = c - 2k z L3 cos(2a^L 3 7 P)C, 7 e- 2K 'r L * cos(2a s P L 2 7 7 )f 

70 fe~^ 2K E + ' 5 E) L 3 4. e -( 2 «P+< 5 p) L 2 _|_ c -2(4L 3 1-K p I 4 ) 



/it = c ~ 2 ^ Ll cos{2a s P L 4 7 7 )£ + ^ e -( 2 ' t p+' 5 p)- L4 + c -2(« E £<3+ K pL2)^ _ 

Subtracting these equations gives the solvability condition 

= {e~ 2K P L ^ cos(2apL 2 + 7) - e" 2 ^ 4 cos(2a s P L 4 + 7)) £ 
+ e -^ u E L 3 C0S ( 2(T « L 3 + 0)C 

+ Q ^ c -(2Kp+«5p)i4 _|_ c -(2Kp+<5p)L 2 _j_ c -(2K^+(5g)L 3 ~j 

+ C ^e _2 ^E i 3+ K P L 4) 7 e -2(KE i 3+K P i2)~j _ 

Since L 2l L 3 ,L 4 are free parameters (for min{L2, L 3 , L4} > L* and within their 
domains) a natural starting point to find solutions are the asymptotic regimes: 

1. (a) L 3 , L 2 > (b) L 3 , L 4 > L 2 ; (c) L 2 , L 4 > L 3 ; 

2. L3 >• L2,Li; 

3. X3 ~ L 2 ~ £4; 

This is at first independent of the relative sizes of Kp " E ; however, these are relevant 
when estimating the constants in the specific meaning of the '3>' symbols. 
Case 1. For subcase (a) the asymptotic solvability condition reads 

= Q-iKpL* cos(2a s P L 4 7 7 )£ + O f e -^+ s ^ L 



For C ^ solutions are close to roots of the cosine-term, which, for up 7^ 0, are 
L 4 = ((2m + 1)tt - j)/2a p , m E N. However, L 4 = iT P /2, £ G N and for generic T P 
the previous fails for any m, I, so that, also for a s P = 0, typically there is no solution 
to the asymptotic equation. In the resonant case Tp = ((2m* 7 l)n — -f)/£*crp, and 
with £ = m£* in L4, the full solvability condition reads 

= e" 2 ^ 2 cos(2apL 2 7 7 )£ + e' 2 ^ 3 cos(2a^L 3 7 /3)( 

j^q ^ e -(2rep+5p)L 4 _|_ e -(2Kp+5p)L 2 _|_ e -(2n E +8 E )L 3 ^j 

This can be solved under the condition L 3 ,L 2 3> L4 if the cosine term with the 
continuous L 3 is leading order. If the difference in l 3>' is not too big, this is possible if 
C 7^ 0, Kp 7 <5 P /2 > Kg, and requires £ = m 2 ^», unless k s p > or L 2 3> L 3 . Looking 
at the remainder term, L 3 3> L 2 is possible if Kg 7 tp > Kg. In all these cases, 
if (Tg 7^ and if L 2 ,L 4 are sufficiently large, the implicit function theorem yields a 
sequence of solutions in L 3 . 

Subcase (b) is analogous, but in subcase (c) L 3 is a continuous parameter so 
that, if £ 7^ and <7g 7^ 0, there exists an infinite sequence of robust solutions to 
the asymptotic equation without the resonance assumption and the constraint on the 
spectral gaps. 
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Case 2. The asymptotic solvability condition reads 

= (c^ 2 ^ 2 cos(2cr^L 2 + 7) - c _2 ^ i4 cos(2a s P L i + 7)) £ 

+ (c-( 2k p+ s p) l 4 + e -(2Kp+«p)ia^ _ 

If £ 7^ and in the non-resonant case, the sum of the cosine-term vanishes if and only if 
L2 = L4. This means that the corresponding orbit revolves about P the same number 
of times each passing. However, discreteness implies that L2, £4 cannot be adjusted 
to compensate the error terms. Hence, we look at the full solvability condition with 
L2 = L4. Since L3 ;§> L4 = L2 this provides a solvable condition if £ ^ and if the 
cosine term in L3 is leading order. Similar to Case l(a),(b), this can be achieved if 
k e < min{ftp + <5p/2, k s e + k s p }, and the result is as in Case 1(c). Hence, if under 
these conditions a long time is spent near E, then the number of rotations about P 
are locked: they must be the same at each passing. 

Case 3. If k e < k s p the situation is as in Case 1(c), and if k e > K S P , typically as 
in Case 2, and under resonance analogous to Case 1(a). 

5.3.2. EPlt. For simply recurrent solutions in the EPlt case, the loci of param- 
eters are given explicitly in Corollary 15.11 but one is also interested in the location 
of turning points in the parameter curves and folds of solutions. In [3] EPlt cycles 
are studied in K 3 using a not entirely rigorous, but geometrically intuitive approach 
to obtain bifurcation equations for simply recurrent solutions. Note that in M 3 cither 
crp U = (positive Floquet multipliers) or <7p U = n/Tp (negative ones) and without 
loss of generality one-dimensional unstable manifold so that a E = 0. It turns out 
that the present rigorous approach confirms the results of [3] in arbitrary ambient 
dimensions. 

Here we briefly illustrate the location of turning points and the bifurcation set 
for 1-homoclinic orbits. 

1-homoclinic orbits to E . This case yields the bifurcation equations 

H! = c~ 2k p l cos(2cr p L + ft + + Civ) 

H2 = av 2 + c - 2k f l cos(2apL + 72)62 ■ 

where v is continuous with \v\ < e and L = £Tp/2 (Tp the period of P) for t sufficiently 
large counting the number of oscillations about P. Up to an error of order v 2 , we can 
write the equation for fii as 

Ml = c - 2 ^ L (cos(2a^i + ft)G + (cos(2a P L + p x )C,[ - sm(2a P L + pi)tift)v). 

Solving for v and substituting the result into the equation for [12 gives 

/ e 2 ^Vi -cos(2a£Z + ft )& \ 2 _ 2k * L /n s T 
^^H cos(2^/ft)C(-sin^ + ft) Cl ft J +e ^cos(2^ + 72 )6, 

note that the denominator typically never vanishes if ^ or (if3[ ^ because 
L = ITp /2 is constrained to a equi-distancc discrete sequence. 
Solving i9 M1 /i2 = gives the turning points of solution curves 

/ii = := c~ 2K r L cos(2o- p L + ft)Ci- 

Hence, the solution set typically is the union of parabolas with critical points at 
/Lti = /i*. Depending on a P and o~ s P the critical points lie on the discrete evaluation of 
either a monotone or a 'snaking' curve in the fix and ^2 direction, respectively, which 
can generate a spiraling sequence in the parameter plane. 
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l-homoclinic orbits to P. The bifurcation equations become 

m = c- 2K - L cos(2o^L + 71 + + £«) 

^2 = ™ 2 + e _2 ^ L cos(2cr^ J L + /3 2 )C2, 

where L > L* and u, < e are both continuous so that solutions typically come in 
two-dimensional sheets connected by folds or corners. 

Eliminating v as in the -E'-homoclinic case gives (for non-resonant L) 

( e 2 ^ L /i! - cos(2g|L + 7l )& 
/i2 a \cos(2a E L + 7l )£ - sin(2^i + 7l )£ l7 J 

Solving d^ ± fi2 = yields the turning points 

A*i = A** :=e- 2K E i cos(2cr| ; L + 7i)^i, 

which give the location of a fold curve of the solution sheet in the parameter plane. 

For each non-resonant sequence of L- values the parameter locations arc analogous 
to those of l-homoclinic orbits to E. If a s E — a E = there are no resonances and the 
monotone /J,*(L) provides all folds. Otherwise, for each resonant value of L, note that 
/X2 is quadratic in v with critical point at v = 0, which means fix = fi*. Hence, to 
leading order, fold curves are given on the one hand via /z*(L), and on the other hand 
via H2{v) at \i\ = \i*(V) for resonant L. A more detailed description of the solution 
set for n — 3 is given in [4]. 



J + c - 2 ^ i cos(2a£L + /3 2 )C 2 . 
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